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Abstract 

This  report  summarizes  the  accomplishments  under  AFOSR  Award  FA9550- 12- 1-0098. 
The  main  goal  of  this  project  was  to  investigate  by  means  of  nanocalorimetry  both  the 
synthesis  of  diboride-based  coatings  using  reactive  multilayers  and  their  oxidation 
behavior.  At  the  onset  of  the  project,  a  scanning  AC  nanocalorimetry  technique  was 
developed  to  make  nanocalorimetry  measurements  possible  at  the  elevated  temperatures 
required  for  the  project.  This  novel  technique  enables  nanocalorimetry  measurements 
over  a  wide  range  of  temperatures  and  scanning  rates  not  accessible  to  other  calorimetry 
techniques.  The  dynamic  range  of  the  technique  makes  it  ideal  for  analyzing  the 
thermodynamics  and  kinetics  of  phase  transformations,  solid-gas  reactions,  and  solid- 
state  reactions  in  thin-film  samples. 

The  scanning  AC  calorimetry  technique  was  employed  to  investigate  the  synthesis  of 
ZrB2  and  carbon-doped  ZrB2  using  Zr/B  and  Zr/B4C  reactive  multilayer  coatings.  The 
solid-state  reactions  in  these  multilayers  were  shown  to  proceed  in  two  distinct  steps:  an 
inter-diffusion/amorphization  step  and  a  crystallization  step.  Measurements  performed  at 
heating  rates  ranging  from  3,100  to  10,000  K/s  allowed  determination  of  the  kinetic 
parameters  of  the  multilayer  reactions,  yielding  activation  energies  of  0.37  -  0.47  eV  for 
the  amorphization  step,  and  1.85  -  2.4  eV  for  the  crystallization  step.  Density  functional 
theoretical  computations  provide  insight  in  the  amorphization  process  and  confirm  the 
relatively  low  activation  energy  associated  with  this  process.  The  presence  of  carbon 
increases  boron  mobility  and  lowers  the  activation  barrier  for  both  steps. 

The  oxidation  resistance  of  diboride-based  coatings  was  evaluated  by  measuring  the  heat 
capacity  of  a  sample  during  the  oxidation  reaction,  a  technique  first  demonstrated  on 
zirconium  thin  films.  Oxidation  of  Zr  was  shown  to  be  diffusion-limited  with  an  energy 
barrier  of  0.59  ±  0.03  eV,  suggesting  a  grain  boundary  diffusion  mechanism.  Oxidation 
of  ZrB2  thin  films  follows  the  Deal-Grove  oxidation  model,  with  activation  energies  of 
0.71  eY  for  the  interfacial  reaction,  and  2.35  eV  for  the  diffusion  step.  At  temperatures 
below  900  K,  oxidation  of  ZrB2  follows  a  logarithmic  rate  law.  In  contrast  to  pure  ZrB2 
films,  carbon-doped  ZrB2  films  experienced  gross  delamination  and  spalling  during 
oxidation  as  a  result  of  gas  formation  during  the  oxidation  process.  Doping  of  ZrB2 
coatings  with  Nb  results  in  a  lowering  of  the  overall  oxidation  resistance  of  the  coatings 
at  temperatures  below  1000  K,  but  the  temperature-dependence  of  the  diffusion  rate 
constant  suggests  that  Nb  doping  may  be  beneficial  at  more  elevated  temperatures. 

Tensile  experiments  were  performed  on  900  nm  sputtered  films  of  un-doped  ZrB2  at 
temperatures  ranging  from  ambient  to  1016  K.  The  ZrB2  films  had  an  ambient- 
temperature  tensile  strength  in  excess  of  1.2  GPa,  far  exceeding  the  tensile  strength  of 
bulk  ZrB2,  but  their  stiffness  was  lower  than  for  bulk  materials  and  decreased 
significantly  with  increasing  temperature.  The  fracture  toughness  of  the  ZrB2  films  was 
2.57  +  0.01  MPaVm,  which  is  similar  to  values  reported  for  bulk  ZrB2. 
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I.  Introduction 

This  report  summarizes  the  results  of  the  research  project  performed  for  the  Air  Force 
Office  of  Scientific  Research  (AFOSR)  under  Contract  FA9550-12-1-0098,  “Fabrication 
and  characterization  of  novel  refractory  coatings  using  combinatorial  nanocalorimetry”. 
The  main  objectives  of  the  project  were  to  investigate  by  means  of  nanocalorimetry  both 
the  synthesis  of  diboride-based  coatings  using  reactive  multilayers  and  their  oxidation 
behavior. 

Ceramics  based  on  zirconium  diboride  have  a  unique  combination  of  very  high  melting 
point,  high  hardness,  great  resistance  to  erosion  and  oxidation,  and  excellent  thermal 
conductivity1.  This  combination  of  properties  makes  them  excellent  candidates  for  use  as 
thermal  protection  materials  in  sharp  leading  edges  of  hypersonic  flight  vehicles  or  in 
nose  cones,  thus  enabling  more  agile  vehicles  that  would  open  up  a  greater  range  of 
hypersonic  flight  paths2.  Bulk  ZrB2-based  ceramics  are  usually  fabricated  from  powders 
that  are  compacted  by  a  range  of  techniques,  including  hot  pressing,  pressure-less 
sintering,  reactive  densification,  and  spark  plasma  sintering.  The  properties  of  bulk 
ceramics  fabricated  using  one  of  these  routes  depend  on  the  processing  conditions  and  on 
the  degree  of  densification  that  is  attained1'2.  Reliable  fabrication  of  bulk  ultra-high- 
temperature  ceramics  is  difficult  and  time  consuming.  By  contrast,  growing  dense 
coatings  of  the  same  materials  by  means  of  sputter  deposition  is  relatively  straightforward. 
Thus  coatings  can  be  used  to  explore  the  properties  of  a  wide  range  of  materials  systems 
in  much  less  time  than  bulk  materials  and  can  provide  valuable  guidance  in  the  design  of 
bulk  materials. 

While  use  of  vacuum  deposition  techniques  greatly  facilitates  the  fabrication  of  ceramic 
samples,  the  small  mass  of  the  samples  imposes  severe  demands  on  the  tools  used  to 
characterize  them.  We  have  developed  a  novel  nanocalorimetry  technique  that  enables 
rapid  screening  and  characterization  of  very  small  quantities  of  materials3'7.  The 
technique  relies  on  a  micromachined  device  with  an  array  of  thermal  sensors  that  can  be 
operated  at  high  temperature.  Calorimetry  can  be  used  to  investigate  the  thermodynamics 
and  kinetics  of  many  phenomena,  ranging  from  phase  transformations  and  reactions  in 
multilayered  films,  to  grain  growth  and  evaluation  of  radiation  damage8'15. 

The  specific  research  goals  of  the  project  include  the  following:  (1)  Fabricate  diboride- 
based  ultra-high-temperature  ceramic  coatings  by  means  of  reactive  multilayers;  (2) 
deposit  coatings  with  in-plane  composition  gradients  by  simultaneous  sputtering  from 
several  different  magnetrons  to  investigate  the  effects  of  additives  on  diboride  oxidation 
behavior  over  a  range  of  compositions;  (3)  Use  scanning  and  isothermal  nanocalorimetry 
techniques  to  investigate  the  oxidation  behavior  of  diboride  coatings  over  a  range  of 
compositions  and  temperatures;  (5)  Characterize  the  mechanical  behavior  of  diboride 
coatings  both  at  room  temperature  and  at  elevated  temperature. 

The  project  report  is  organized  as  follows:  In  Section  III,  we  provide  a  brief  description 
of  the  device  that  was  used  to  perform  the  nanocalorimetry  measurements  and  we  discuss 
the  theory  for  scanning  AC  nanocalorimetry  that  was  developed  as  part  of  this  project.  In 
Section  IV,  we  describe  the  tensile  tester  that  was  used  to  determine  the  mechanical 
behavior  of  ZrB2  thin  films  as  a  function  of  temperature.  In  the  Section  V,  we  present  the 
results  obtained  for  several  materials  systems.  The  conclusions  are  presented  in 
Section  VI.  The  Appendices  contain  the  descriptions  of  several  measurement  procedures 
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that  were  developed  as  part  of  this  project,  as  well  as  an  example  how  the  scanning  AC 
nanocalorimetry  technique  can  be  applied  to  the  broader  field  of  materials  science. 

II.  Dissemination  of  results 

The  results  of  this  research  project  have  led  to  12  publications  in  archival  journals  with 
several  more  in  preparation,  to  10  invited  or  keynote  presentations,  and  several 
contributed  conference  presentations.  The  publications  are  listed  in  Section  VIII. 

III.  Scanning  nanocalorimetry 

In  this  section  we  provide  a  brief  description  of  the  sensors  (Fig.  1)  that  were  used  to 
perform  the  nanocalorimetry  measurements.  We  also  present  the  theory  underlying 
scanning  AC  nanocalorimetry,  a  technique  developed  as  part  of  this  program  that  has 
made  it  possible  to  perform  accurate  calorimetry  measurements  at  elevated  temperatures 
with  a  dynamic  range  spanning  many  decades  of  heating  rate. 

A.  Sensors 

The  nanocalorimeter  device  consists  of  a  silicon  substrate  with  an  array  of  micro- 
machined  thermal  sensors.  The  device  is  fabricated  using  conventional  Si-based 
microfabrication  steps  as  illustrated  in  Fig.  2.  The  sensors  are  arranged  in  a  5  x  5  array  to 
facilitate  combinatorial  sample  preparation  and  analysis  (Fig.  1).  When  a  thin-film 
sample  with  an  in-plane  composition  gradient  is  deposited  on  the  substrate,  the  film  is 
essentially  discretized  at  each  sensor,  allowing  the  simultaneous  creation  of  25  samples 

with  unique  composition.  Each  sensor 
consists  of  a  thin-film  thermistor  sand¬ 
wiched  between  two  electrically 
insulating  silicon  nitride  layers  that  form 
a  membrane  supported  by  the  substrate 
(Fig.  3(a)).  The  thermistor  serves  both  as 
a  heating  element  and  a  thermometer.  It 
is  made  of  tungsten  because  of  its  large 
temperature  coefficient  of  resistance  and 
its  small  resistivity,  both  of  which  are 
beneficial  to  measurement  sensitivity. 
The  high  melting  temperature  of  tungsten 
also  results  in  excellent  thermal  stability 
of  the  thermistor.  The  electrical  leads  and 
Fig.  l.  Photograph  of  the  PnSC  device  contact  pads  on  the  substrate  are  made  of 

copper  to  reduce  the  resistance  of  the 
signal  lines  on  the  substrate  and  to  facilitate  contact  to  the  PnSC  device. 

Samples  to  be  measured  are  limited  to  the  thermistor  area  of  each  sensor  through  use  of  a 
shadow  mask.  The  membrane  design  of  the  sensor  thermally  insulates  the  sample  from  its 
surroundings  and  ensures  that  the  thermal  mass  of  the  sensor,  i.e.,  the  addendum,  is  very 
small.  Referring  to  the  schematics  in  Fig.  3,  the  straight  lines  down  the  center  of  the 
membrane  is  the  heating  element;  the  metal  lines  connected  to  the  heater  are  the  voltage 
probes,  and  the  portion  of  the  heating  element  between  the  voltage  probes  is  the 
thermistor. 
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Fig.  2.  (a)  The  fabrication  process  starts  with  a  (100)  Si  wafer  passivated  with  300  nm  low-stress 
LPCVD  nitiride.  (b)  Photoresist  is  spin-coated  and  patterend  on  the  front  side  of  the  wafer,  (c)  Layers 
of  5  nm  titanium  and  100  nm  tungsten  are  sputter  deposited  and  patterned  using  a  lift-off  procedure. 
Copper  is  then  deposited  and  patterned  using  the  same  procedure  (not  shown  in  the  figure)  to  connect 
the  heater  to  the  contact  pads,  (d)  After  metal  patterning,  30  nm  Si3N4  is  coated  on  the  front  side  to 
protect  the  metal  layers  from  oxidation,  (e)  To  make  freestanding  membranes,  photoresist  is  patterned 
on  the  backside  and  (f)  Si  is  removed  using  a  SPTS  Rapier  DRIE  system. 


A  current  passing  through  the  thermistor  heats  the  sample  and  the  calorimetric  cell.  The 
power  dissipated  in  the  thermistor  is  determined  experimentally  from  the  current  supplied 
to  the  thermistor  and  the  potential  drop  between  the  voltage  probes.  The  local 
temperature  change  is  determined  from  a  four-point  thermistor  resistance  measurement 
that  has  been  calibrated  to  temperature.  The  calorimetric  signal  is  then  calculated  from 
the  power  dissipated  in  the  thermistor  and  the  heating  rate  of  the  sample. 

Not  all  power  dissipated  in  the  thermistor  is  used  to  heat  the  sample  and  addendum;  some 
of  it  is  lost  to  the  environment.  At  moderate  temperatures,  heat  transfer  from  the 
thermistor  to  the  membrane  dominates  this  heat  loss.  As  the  temperature  increases, 
radiation  from  the  thermistor  and  the  membrane  becomes  important  and  eventually 
dominates  the  heat  loss.  Measurements  are  generally  performed  in  vacuum  to  eliminate 
convection  losses  and  to  provide  a  chemically  inert  testing  environment.  Conduction  and 
radiation  heat  losses  can  be  accounted  for  by  modeling  or  through  the  use  of  a  reference 
measurement  scheme. 

Alternatively,  measurements  can  be  performed  using  scanning  AC  calorimetry  method, 
where  measurement  results  are  insensitive  to  heat  loss.  In  this  project,  a  combination  of 
scanning  DC  and  AC  calorimetry  methods  were  employed  to  study  the  materials  systems. 

B.  Scanning  AC  nanocalorimetry 

AC  techniques  have  long  been  used  in  bulk  calorimetry  measurements  where  the 
underlying  temperature  is  ramped  up  very  slowly16'19,  the  application  of  AC  techniques  in 
fast  scanning  nanocalorimetry  is  new.  The  purpose  of  this  section  is  to  provide  the 
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Fig.  3.  Layout  of  the  nanocalorimeter  cell:  (a)  cross-section  schematic,  and  (b)  plan-view  schematic. 
Heater  line  width  is  0.8  mm  and  voltage  probe  line  width  is  0.1  mm. 


underlying  theory  and  the  conditions  under  which  valid  AC  measurements  can  be 
performed  using  fast  scanning  nanocalorimetry. 

Bl.  Analysis  of  thermistor-based  sensors 

Consider  a  nanocalorimetry  sensor  that  consists  of  a  heating  element  with  a  temperature- 
dependent  resistance.  The  sample  of  interest  is  placed  on  the  sensor  and  is  heated  by 
sending  a  current  through  the  heating  element.  The  temperature  of  the  sensor  is 
determined  by  monitoring  the  resistance  of  the  heater  using  a  four-point  measurement 
scheme  (Fig.  4),  while  the  power  supplied  to  the  sample  is  determined  from  Joule’s  law. 
The  following  analysis  assumes  that  the  temperature  is  uniform  across  the  heating 
element  and  that  the  temperature  difference  between  the  heating  element  and  the  sample 
is  negligible.  The  first  assumption  depends  on  the  planar  geometry  of  the  sensor.  See 
reference  [20]  for  a  more  detailed  discussion  of  temperature  uniformity.  The  second 
assumption  is  valid  as  long  as  the  thermal  diffusion  length  is  significantly  larger  than  the 
sensor  thickness  at  the  time  scale  of  the  measurements,  which  is  indeed  the  case  for  the 
thermistor-based  sensor  if  the  AC  frequency  is  not  too  large. 

In  a  typical  AC  nanocalorimetry  measurement,  a  current  with  both  DC  and  AC 
components, 


V 


Sample 

Heater 


Fig.  4.  Schematic  of  a  four-point  measurement  on  a  heating  element  with  sample. 
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/  =  70  +  i  cos(<x>t  —  <p0),  (3.1) 

is  passed  through  the  resistive  heater.  The  AC  component  causes  the  temperature  of  the 
heater  to  oscillate  around  an  average  temperature  T0  that  changes  with  time.  Hence  the 
temperature  response  of  the  heater  can  be  described  as 

7  =  70(0 +  0(0,  (3.2) 

where  T0(t)  is  the  average  temperature  over  one  oscillation  period  and  0(t )  represents  the 
oscillating  part  of  the  temperature.  Since  the  resistance  of  the  heating  element  is 

temperature  dependent,  the  temperature  oscillations  cause  the  resistance  of  the  heating 

21 

element  R  to  vary  as 

R=R0(1+A06 ),  (3.3) 

where  Ro  and  Aq  are  the  resistance  and  the  thermal  coefficient  of  resistance  (TCR)  of  the 
heating  element  at  To,  respectively.  The  power  dissipated  by  the  heater  follows  from 
Joule’ s  law 

P  =  I2R  =  [I0  +  i  cos  (cot  -  (p0)]2R0(l  +  A06).  (3.4) 

Expanding  the  rate  of  heat  loss  from  heater  and  sample  in  a  Taylor  series  around  To  yields 

to  first  order 

L(T)  =  L0(T0)  +  a06,  (3.5) 

where  Lo(7o)  and  a()  are,  respectively,  the  rate  of  heat  loss  and  its  temperature  derivative 
at  To.  In  the  absence  of  a  phase  transition  or  a  chemical  reaction,  the  energy  balance  takes 
the  form, 

P  —  CT  +  L(T)  =  Cp(To  +  0^  +  Lq  +  (XqO  .  (3.6) 

where  C  is  the  total  heat  capacity  of  sample  and  addendum  (Fig.  3).  To  ensure  that  Ro 

remains  constant  on  the  time  scale  of  the  oscillations,  it  is  sufficient  that 

A0T0  —  «  1.  (3.7) 

CO 

Because  both  A  and  0  are  small,  it  further  follows  that 

A06>  «  1.  (3.8) 

A  combination  of  Eqs.  (3.4)  through  (3.8)  then  yields 

Oo  +  if2)£0  +  2 10iR0  cos  (cot  -  <p0)  +  \i2Ro  cos2  (<ut  -  <p0)  =  CP(T0  +  G)  + 

Lq  +  a00.  (3.9) 

If  we  take  the  time  average  of  Eq.  (3.9)  over  one  oscillation  period,  all  oscillating  terms 
vanish,  yielding  the  following  equation  in  To, 

QZ+li2)R0  =  CpT0+L0.  (3.10) 

Eq.  (3.10)  provides  the  relationship  between  the  DC  component  of  the  applied  power,  the 
heat  loss,  and  the  ramp  rate  for  the  average  temperature  To.  Subtracting  Eq.  (3.10)  from 
Eq.  (3.9)  leads  to  a  first-order  ordinary  differential  equation  in  6, 

2 10iR0  cos (rnt  —  <p0)  +^i2R0  cos2 (cut  —  <p0)  =  Cv6  +  a00.  (3.11) 
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This  equation  describes  the  oscillating  component  of  the  temperature.  The  temperature  To 
enters  into  the  equation  only  via  the  temperature-dependence  of  Ro  and  a(),  which  can  be 
regarded  as  constant  on  the  time  scale  of  the  oscillations  by  Eqs.  (3.5)  and  (3.7).  The 
steady-state  solution  to  Eq.  (3.11)  is  given  by 


6  =  9±  cos (wt  -  cp0  -  (Pi)  +  02  cos(2 cot  -  2 cp0  -  cp2), 


0i  = 


9 ,  = 


HIqRq  sin 
Coo 

i2Rp  sin  cp2 
4  Coo  ’ 


tan<p!  =  -tan<p2 


Coo 

a0' 


(3.12) 

(3.13) 

(3.14) 

(3.15) 


While  the  current  passing  through  the  heating  element  has  an  angular  frequency  co,  the 
temperature  oscillation  has  both  co  and  2 co  components  because  the  power  scales  with  the 
square  of  the  current.  The  amplitudes  of  both  components  scale  inversely  with  the  heat 
capacity  C;  the  phase  angles  of  the  temperature  oscillations  with  respect  to  the  applied 
AC  current  are  increasing  functions  of  the  ratio  Cco/cXq.  Once  the  temperature  response  of 
the  sensor  is  known,  the  voltage  response  of  the  thermistor  is  easily  derived.  Substituting 
the  solution  for  0  given  by  Eq.  (3.12)  into  Eq.  (3.3)  provides  an  expression  for  the 
resistance  of  the  heating  element, 

R  =  R0[  1  +  A061  cos  (cot  -  (Po  ~  <Pi)  +  A092  cos(2  cot  -  2  cp0  -  cp2)  ]  •  (3.16) 

The  voltage  drop  V across  the  heating  element  is  then  given  by  Ohm’s  law, 

V  =  [I0  +  i  cos (cot  —  cp0)]R0[l  +  A091  cos(o>t  —  <p0  —  cp1 )  +  A092  cos(2 cot  —  2 cp0  — 
cp2)\.  (3.17) 


In  addition  to  a  DC  component,  the  voltage  response  of  the  sensor  has  co,  2  co  and  3  co 
harmonic  components.  It  is  useful  to  deconvolve  the  voltage  response  into  the  four 


frequency  components,  Vo,  V»,  V2a,  and  V?®.  The  DC  component  is  given  by 

Vo  =  IoRo  T  ~  cos  <Pi  ~  I0R0,  (3.18) 

while  the  harmonic  components  are 

cos  (ot  -1-  sin  cot,  (3.19) 

=  iR0  cos  (p0  +  I0R0A091  cos(<p0  +  <Pi)  +  \  iR0^ 0^2  cos (<p0  +  <p2\  (3.20) 

Ya>  =  iR0  sin <p0  +  I0R0A09 1  sin (cp0  +  <p±) +\ iRoA0G2  sin(<p0  +  <p2),  (3.21) 

V2 00  -  X200  cos  2<x)t  +  Y2o)  sin  2 cot,  (3.22) 

X2w  =  \ iff o^tA  cos(2<p0  +  <pt)  +  IR0A092  cos(2 <p0  +  <p2),  (3.23) 

y2oi  =  \ iff o^offi  sin(2<p0  +  <Pi)  +  IR0A0d2  sin(2 <p0  +  cp2),  (3.24) 

^3a>  =  X3co  cos  3 cot  +  YlCO  sin  3  cot,  (3.25) 

X3to  =  \ iRoA0e2  cos(3<p0  +  (pz),  (3.26) 
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(3.27) 


*30.  =  \iR~ o^oe2  sin(3^0  +  <P2)- 

Once  the  resistance  of  the  heating  element  is  calibrated  to  temperature,  To  can  be 
determined  directly  from  Eq.  (3.18).  According  to  Eq.  (3.8),  the  terms  containing  6\  and 
62  in  Eq.  (3.20)  and  (3.21)  are  much  smaller  than  terms  in  iRo,  so  that  they  do  not 
generally  provide  accurate  information  on  the  heat  capacity.  It  is  possible,  however,  to 
determine  the  heat  capacity  from  either  Eqs.  (3.22)  through  (3.24),  or  Eqs.  (3.25)  through 
(3.27).  Considering  that  6\  is  usually  more  than  an  order  of  magnitude  larger  than  62 ,  the 
2co  component  provides  a  stronger  signal  than  the  3©  component.  Substituting  the 
expressions  for  0\  and  G2  into  Eqs.  (3.23)  and  (3.24)  yields 

x2a>  =  '-^±[sm<p1cosO.cp0  +  <p1)+-sm<p2cos(2cp0  +  cp2)],  (3.28) 

LpCO  4 

Y2o,  =  ‘  [sin  cp1  sin(2(p0  +  cp2)  +  ^sin  <p2  sin(2<p0  +  (p2)].  (3.29) 

LpUJ  4 

Equations  (3.15),  (3.28)  and  (3.29)  can  be  solved  for  (p\,  q>2,  C  and  a.  The  phase  angle  (p\ 
is  given  by  the  following  equation,22 

10  tan3  (p1  +  9 n  tan2  <p±  +  4  tan  +  3n  =  0,  (3.30) 

where 


_  X2a)sm2(p0-Y2(l)cos2(p0  (3  31) 

X20)  cos  2cp0 +Y2a)  sin  2<p0' 

Equation  (3.30)  has  only  one  real  root,  so  that  tan  (p\  is  uniquely  determined.  The  heat 
capacity  of  the  sample  and  addendum  is  then  given  by22 


i2lQRQ?iQ  (25  tan2  (p±+9)  tan2  cp1 
u>\V2o)\  J  (16tan2  <p1+4)(l+tan2  cp^’ 


(3.32) 


Under  a  broad  range  of  experimental  conditions,  Cco  is  much  larger  than  a.  According  to 
Eq.  (3.15),  the  phase  angles  (p\  and  qh  then  approach  a  value  of  tt/2  and  Eq.  (3.32) 
reduces  to  the  following  simple  form, 


5i2I0R^0 

4w|V2&)| 


(3.33) 


If,  on  the  other  hand,  Cco  is  much  smaller  than  a,  evaluation  of  Eq.  (3.32)  is  prone  to 
large  relative  errors  and  C  is  difficult  to  determine  accurately.  Instead  of  performing  a 
harmonic  analysis  of  the  voltage  signal,  it  is  in  principle  also  possible  to  calculate  C  and 
a  from  an  analysis  of  the  resistance  of  the  thermistor.  The  parameters  of  interest  can  then 
be  calculated  from  using  Eq.  (3.16).  At  first  glance,  this  approach  seems  simpler  than  the 
method  outlined  above,  but  this  approach  works  poorly  when  70  and  i  have  similar 
magnitudes  because  of  the  poor  signal-to-noise  ratio  when  the  instantaneous  applied 
current  approaches  zero. 

B2.  Selection  of  experimental  parameters 

Equations  (3.7)  and  (3.8)  represent  the  conditions  under  which  a  scanning  AC 
calorimetry  measurement  can  be  analyzed  using  Eq.  (3.32).  When  designing  an 
experiment,  parameters  such  as  the  DC  and  AC  current  components  and  the  AC 
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frequency  need  to  be  selected  to  satisfy  these  conditions.  Substituting  Eqs.  (3.12)-(3.14) 
into  Eq.  (3.8)  leads  to  the  following  equivalent  condition20, 


xe  0  +  ij)  =  ct « l. 

CpO)  V  U  4  J  1 


(3.34) 


In  addition  to  these  conditions,  there  are  two  experimental  considerations:  1)  the 
amplitude  of  the  2 co  voltage  signal  should  be  at  least  20  times  greater  than  the  resolution 
of  the  voltage  measurements;  2)  the  angular  frequency  co  must  be  large  enough  that  the 
mean  temperature  rise  in  a  given  cycle  does  not  exceed  the  desired  temperature 
resolution,  typically  about  2  K1.  These  considerations  lead  to  the  following  conditions20, 


V2U  _  si2i0r^a0  _  20 

?  4  CpO)% 


(3.35) 


and 


T0-=C3<  2K,  (3.36) 

CO 

where  £  represents  the  resolution  of  voltage  measurements.  Table  I  lists  typical  values  of 
experimental  parameters  for  a  scanning  AC  nanocalorimetry  measurement.  Note  that  for 
most  thermistor  materials,  Eq.  (3.7)  is  automatically  satisfied  if  Eq.  (3.36)  is,  and  only 
conditions  (3.34)-(3.36)  need  to  be  considered  when  selecting  the  experimental 
parameters. 

If  the  measurement  is  conducted  in  an  environment  with  high  heat  loss,  the  sensor 
undergoes  a  quasi  steady-state  scan  and  the  value  of  T0  is  determined  by  the 
instantaneous  values  of  7o  and  i.  The  nominal  scanning  rate,  t0,  is  then  controlled  by  the 
ramp  rates  of  7o  and  i  only,  and  can  be  selected  independently.  Once  the  scanning  rate  is 
set,  the  value  of  the  angular  frequency  is  selected  such  that  condition  (3.36)  is  satisfied. 
Allowable  combinations  of  7o  and  i  are  then  determined  from  conditions  (3.34)  and 
(3.35). 

Fig.  5  illustrates  a  simple  graphical  method  for  selecting  the  current  parameters.  The 
shaded  area  in  the  figure  represents  the  locus  of  permissible  combinations  of  7o  and  i  for 
typical  experimental  parameters  (Table  I).  The  boundary  of  the  shaded  area  depends 
slightly  on  temperature  -  as  the  temperature  increases  from  350  K  to  870  K,  the 
permissible  locus  increases  in  size  as  indicated  in  the  figure.  The  two  red  lines  are 
contours  of  constant  average  power.  If  the  sensor  is  operated  under  quasi  steady-state 
conditions,  they  also  represent  isotherms.  Connecting  two  points  within  the  permissible 
region,  one  on  each  of  the  isotherms,  with  a  straight  line  provides  a  locus  of  possible 
current  combinations  that  leads  to  a  temperature  rise  from  350  K  to  870  K.  A  current 
program  that  follows  this  path  through  parameter  space  provides  a  quasi  steady-state 


Table  I.  Typical  experimental  parameters 


C  (J/K) 

^RT  (£2) 

To  (K) 

A  (1/K) 

t0  (K/s) 

100 

a  (W/K) 

2  -  8><10‘6 

5-  10 

370-1400 

1.3  -  1.6xl0‘3 

101-  103 

10"5 

lO^.lO’3 
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Fig.  5.  Graphical  solution  to  inequality  set  (3.34)  and  (3.35).  Condition  (3.34)  is  independent  of 
temperature,  while  condition  (3.35)  shifts  with  increasing  temperature  as  marked  in  the  figure.  The 
curves  of  constant  input  power  correspond  to  the  350  K  and  870  K  isotherms. 


sweep  in  temperature  between  the  two  extremes.  Thus,  the  current  profile  supplied  to  the 
sensor  depends  not  only  on  the  scanning  rate,  but  also  on  the  particular  path  chosen 
within  the  permissible  region. 

For  a  measurement  at  a  controlled  heating  rate,  it  is  necessary  to  select  a  current  profile 
that  maintains  a  monotonically  increasing  power.  Conversely,  an  experiment  at  a 
controlled  cooling  rate  can  be  conducted  by  choosing  a  path  inside  the  permissible  region 
that  maintains  a  monotonically  decreasing  power.  In  a  vacuum  environment,  where  the 
heat  loss  is  much  smaller,  the  quasi  steady-state  assumption  is  only  valid  at  high 
temperatures  (>  600  K).  Control  of  the  scanning  rate  following  the  above  mentioned 
parameter  scan  is  then  only  possible  at  high  temperatures.  At  low  temperatures,  the 
heating  rate  is  determined  by  the  instantaneous  input  power,  which  needs  to  be  chosen  so 
that  Eq.  (3.36)  is  maintained. 

IV.  Micro-tensile  tester 

We  have  built  a  mechanical  tester  for  use  inside  the  SEM  to  investigate  the  mechanical 
behavior  of  thin  films  at  elevated  temperatures.  The  displacement  resolution  of  the  tester 
is  10  nm  with  a  250  pm  stroke;  the  load  resolution  is  9.7  pN  (Fig.  6  (a)).  The  temperature 
of  the  thin-film  samples  can  be  controlled  by  custom  designed  micromachined  heaters 
integrated  with  the  tensile  specimen  (Fig.  6  (b)).  Each  specimen  consists  of  a  tungsten 
heating  element  that  also  serves  as  a  temperature  sensor  and  a  freestanding  thin  film  for 
mechanical  measurement.  Finite  element  modeling  was  performed  to  ensure  a  uniform 
temperature  distribution  in  the  freestanding  section  of  the  tensile  sample  and  to  estimate 
the  temperature  during  Joule  heating.  This  system  was  used  to  characterize  the 
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Fig.  6  (a)  Microheater  attached  to  the  in-situ  tensile  tester  (b)  SEM  image  showing  the  tungsten 
heater  and  the  ZrB2  freestanding  film 

mechanical  behavior  of  thin  ZrB2  films  as  a  function  of  temperature.  The  results  are 
discussed  in  Section  V.C. 


V.  Results 

A.  Synthesis  of  ZrB2-based  coatings  using  reactive  multilayers 
Al.  Nanocalorimetry  and  TEM  studies  of  Zr/B  multilayers 

ZrB2-based  alloys  are  generally  synthesized  via  reduction  processes  ,  chemical  routes  '  , 
or  reactive  processes26,27.  Solid-state  alloying  of  Zr  and  B  is  the  intrinsic  process 
underlying  nearly  all  fabrication  routes.  To  reach  a  comprehensive  understanding  of  the 
solid-state  alloying  process  of  Zr-B,  we  have  investigated  the  energetics,  diffusion 
kinetics,  as  well  as  the  structural  evolution  of  Zr/B  multilayers  using  nanocalorimetry  and 
transmission  electron  microscopy.  Zr/B  multilayer  samples  with  bilayer  periods  of  4.1 
and  8.3  nm  were  sputter  deposited  onto  nanocalorimetry  sensors.  All  multilayer  samples 
in  this  study  had  a  total  thickness  of  66  nm  and  the  thicknesses  of  the  individual  layers 
comprising  the  multilayers  were  chosen  to  ensure  formation  of  stoichiometric  ZrB2  upon 
completion  of  the  reaction. 

Figs.  7  (a,b)  illustrate  the  enthalpy  rate  per  unit  area  determined  from  nanocalorimetry 
measurements.  The  figures  show  the  effect  of  both  heating  rate  and  bilayer  thickness.  The 
scans  reveal  two  separate  exothermic  reactions:  a  reaction  with  a  large  energy  release  and 
a  reaction  with  a  much  smaller  energy  release  at  a  slightly  higher  temperature.  Fig.  7(a) 
shows  that  both  peaks  shift  to  higher  temperatures  as  the  samples  are  heated  faster.  Fig. 
7(b)  clearly  demonstrates  that  the  two  peaks  become  more  distinct  as  the  bilayer  period  is 
reduced.  The  structural  evolution  of  the  multilayer  samples  in  the  various  reaction  stages 
was  investigated  with  cross-section  TEM  as  shown  in  Fig.  7(c). 

The  micrographs  were  obtained  from  the  samples  heated  to  various  temperatures  as 
indicated  in  Fig.  7(b)  before  being  quenched.  The  bright-field  (left)  and  diffraction 
images  (right)  in  Fig.  7(c)  clearly  demonstrate  the  evolution  of  the  samples. 

As-deposited  samples  have  a  structure  that  consists  of  alternating  layers  of  amorphous  B 
and  nano-crystalline  Zr  with  a  (1010)  texture  (Powder  Diffraction  File  #050665).  The  Zr 
and  B  layers  in  the  samples  heated  to  723  and  923K  have  inter-diffused  and  have  started 
to  form  an  amorphous  Zr/B  alloy.  At  723  K,  the  layered  structure  of  the  sample  is  still 
present,  but  the  Zr  diffraction  pattern  is  more  diffuse,  indicative  of  the  onset  of 
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amorphization.  At  923  K  the  layered  structure  has  mostly  disappeared  and  the  diffraction 
pattern  is  now  a  diffuse  halo,  typical  for  an  amorphous  material.  A  few  weak  diffraction 
spots  can  still  be  discerned,  indicating  that  the  amorphization  process  is  not  yet  complete 
at  this  temperature. 

The  sample  heated  to  1,273  K,  above  the  temperature  at  which  the  second  reaction  is 
complete,  has  a  completely  different  structure:  it  consists  of  crystalline  ZrB2  with  a 
(0001)  texture  (Powder  Diffraction  File  #340423).  Thus,  the  second  peak  in  the  enthalpy 
rate  curve  corresponds  to  a  crystallization  reaction.  The  crystalline  ZrB2  has  a  (0001) 
texture,  indicating  that  the  crystalline  phase  nucleated  at  the  top  or  bottom  interfaces  of 
the  multilayer  coating,  or  possibly  at  the  original  locations  of  the  Zr/B  interfaces. 

The  enthalpy  rate  curves  in  Fig.  7(a,b)  can  be  used  to  determine  the  kinetic  parameters  of 
the  two  reactions.  We  first  derive  an  expression  for  the  diffusion  activation  energy 
corresponding  to  the  low-temperature  peak.  In  the  model,  we  assume  that  crystalline  Zr  is 

(a)  (c) 
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Fig.  7.  (a)  Enthalpy  rate  per  unit  area  vs.  temperature  for  Zr/B  multilayer  samples  with  a  bilayer  period  of  8.3nm. 
(b)  Enthalpy  rate  per  unit  area  for  multilayer  samples  with  two  different  bilayer  periods  obtained  at  a  heating  rate 
of  10,000  K/s.  The  temperatures  from  which  samples  were  quenched  for  TEM  observations  have  been  marked  with 
arrows,  (c)  TEM  observations  on  an  as-deposited  sample  (a-c),  a  sample  heated  to  723  K  (d-f),  a  sample  heated  to 
923  K  (g-i),  and  a  fully  reacted  sample  (j-1):  (right)  Bright-Held  images  (right)  Selected  area  diffraction  patterns  at 
each  temperature.  All  samples  had  a  bilayer  period  of  8.3  nm  and  were  heated  to  the  relevant  temperature  at  a 
nominal  rate  of  10,000  K/s,  before  being  quenched  to  room  temperature. 
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Fig.  8.  (a)  Schematic  free  energy  diagram  of  Zr-B  system.  Crystalline  Zr  is  in  meta-stable 
equilibrium  with  an  amorphous  phase  with  B  concentration  Ci.  (b)  Concentration  profiles  after 
inter-diffusion. 


in  metastable  equilibrium  with  an  amorphous  B/Zr  phase  with  boron  concentration  C\. 
The  solid  solubility  of  B  in  c-Zr  is  taken  to  be  negligible  (Fig.  8(a)).  After  inter-diffusion, 
the  phase  boundary  has  moved  a  distance  x,  and  the  B  concentration  in  the  amorphous 
phase  has  become  non-uniform  (Fig.  8(b)).  We  only  consider  the  initial  stage  of  the 
process  where  the  diffusion  distance  is  much  smaller  than  the  bilayer  period  of  the 
sample  and  formulate  the  problem  as  a  moving-boundary  problem  of  the  Stefan-type  . 
The  concentration  profile  is  the  solution  of  the  Fick’s  law  with  the  following  boundary 
conditions, 


C(-°V)  =  1, 

c&t)=cv 


(5.1) 


The  mass  balance  at  the  interface  between  the  amorphous  B  and  the  crystalline  Zr  phases 
requires  that 


__I)  dC 
dt  C1  dx 


(5.2) 


Under  isothermal  conditions,  the  solution  of  this  problem  is  a  variant  of  Neumann’s 
similarity  solution  and  is  given  by 


l-Ci 
l  +  erf  a 


*  f 

l  +  erf 

V  V' 


•.Jot 


(5.3) 


t;  =  2a\H)t , 


(5.4) 


where  D  is  the  inter-diffusivity  and  ms  defined  by 


a 


(l  +  erf  ajexpor 


l-Ci 

4nCi 


(5.5) 
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Taking  the  enthalpy  released  during  inter-diffusion  to  be  proportional  to  the  Zr 
concentration  in  the  amorphous  phase,  it  follows  that 

H  =  —  j(i-c)dx  =  2^-4Dt  ,  (5.6) 

fzr  _oo  tzr 


so  that 


HdH_JaAHX 

dt  {  tZr  J 


(5.7) 


In  this  equation  A H  is  the  total  enthalpy  released  during  inter-diffusion  of  a  single  bilayer 
and  tzx  is  the  thickness  of  the  Zr  layers.  Eq.  (5.7)  can  be  employed  to  determine  the 
diffusivity  from  the  reaction  enthalpy  rate  obtained  in  a  calorimetry  measurement 
performed  under  isothermal  conditions. 

If  measurements  are  not  conducted  under  isothermal  conditions,  the  solution  to  the 
diffusion  problem  (5.3)  is  still  valid,  provided  the  product  Dt  is  replaced  by  a  quantity 
that  depends  on  the  entire  temperature  history  T(t)  of  the  sample, 

T  =  |£>(r(t))dt.  (5.8) 


For  an  experiment  performed  at  constant  heating  rate  /?,  Eq.  (5.6)  then  becomes 


H  =  2 


c/AH 
*Z r\[P 


(5.9) 


but  the  product  HdH/dt  is  still  given  by  Eq.  (5.7)  despite  the  non-isothermal  conditions. 
If  the  diffusivity  D  is  described  by  an  Arrhenius  equation,  it  follows  that 


In 


H 


dH 

dt 


=  2ln 


ccAhJiD^ 


lZr 


kBT ' 


(5.10) 


where  kB  is  the  Boltzmann  constant,  Ed  is  the  diffusion  activation  energy,  and  D0  is  the 
pre-exponential  factor  of  the  diffusion  coefficient.  If  some  intermixing  has  taken  place 
during  deposition  of  the  reactive  multilayer,  Eq.  (5.10)  needs  to  be  replaced  by 


In 


(H  +  H0 


dH 

dt 


=  2ln 


ocAHyfiw^ 


lZr 


A. 

kBT’ 


(5.11) 


where  H„  is  the  enthalpy  associated  with  intermixing.  According  to  Eq.  (5.11),  a  graph  of 
ln((H  +  H0)dH/dt )  versus  1  !kBT  produces  a  straight  line  with  as  slope  the  activation 
energy  Ed. 

Fig.  9(a)  depicts  a  graph  of  In ((H  +  Ha)dH/dt )  versus  I  IksT  for  different  heating  rates. 
All  curves  are  very  nearly  linear,  indicating  that  the  kinetics  of  the  first  reaction  step  is 
consistent  with  a  diffusion-controlled  process,  in  agreement  with  the  TEM  observations. 
The  diffusion  activation  energy  can  be  determined  from  the  slopes  of  the  curves  and  is 
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Fig.  9.  Kinetic  analysis  on  the  reaction  in  Zr/B  multilayer:  (a)  ln((H  +  H0)dH/dt )  versus  1  /kBT  for 
the  temperature  range  of  580-  890  K.  The  slopes  in  this  plot  correspond  to  the  activation  energy  of  the 
diffusion  process,  (b)  Kissinger  plots  for  the  crystallization  peaks  from  the  measurements  on  samples 
with  4.1  and  8.3  nm  bilayer  periods. 


estimated  to  be  0.47  +  0.08  eV.  This  value  is  much  smaller  than  the  1.69-  1.99  eV 
obtained  by  Samsonov  and  colleagues29  for  the  diffusion  of  B  in  Zr.  The  activation 
energy  obtained  by  Samsonov  was  determined  by  measuring  the  thickness  of  the 
crystalline  ZrB2  in  Zr/B  diffusion  couples  exposed  to  temperatures  in  excess  of  1,373  K. 
No  amorphous  Zr/B  phase  was  observed  in  their  experiments.  Thus  formation  of  ZrB2 
under  these  conditions  requires  B  to  diffuse  through  the  crystalline  ZrB2  layer,  a  process 
that  is  evidently  much  more  difficult  than  diffusion  through  the  amorphous  Zr/B  phase. 

The  activation  energy  for  the  crystallization  process  can  be  determined  from  a  Kissinger 
plot  of  the  heating  rate  versus  peak  temperature.  The  results  are  shown  in  Fig.  9(b).  The 
data  for  the  samples  with  the  8.3  nm  bilayer  period  show  a  significant  amount  of  scatter 
because  the  overlap  of  both  peaks  in  the  nanocalorimetry  traces  makes  it  difficult  to 
determine  an  accurate  value  of  the  peak  temperature.  The  results  for  the  multilayers  with 
the  4.1  nm  period  are  better  and  provide  an  estimate  of  the  activation  energy  of  2.4  eV. 
Although  it  is  not  possible  to  get  an  accurate  value  of  the  activation  energy  from  the  data 
points  for  the  8.3  nm  samples,  the  results  are  certainly  consistent  with  this  value. 

A2.  Nanocalorimetry,  XRD,  and  TEM  studies  on  Zr/B4C  multilayers 

Utilization  of  additives  can  improve  the  poor  processability  of  dense  ZrB2  based  alloys. 
Among  the  additives,  carbon  is  well  known  to  be  an  effective  sintering  aid.  It  has  been 
reported  that  presence  of  carbon  promotes  removal  of  surface  oxides  in  ZrB2  and  limits 
grain  growth  of  the  refractory  material,  enhancing  its  sinterability31'34.  However,  the 
beneficial  role  of  carbon  in  the  formation  of  ZrB2-based  alloys  is  still  not  fully 
understood  .  We  have  investigated  the  kinetic  role  of  carbon  in  the  formation  of  ZrB2 
based  alloys.  Zr/B4C  multilayer  samples  were  prepared  for  nanocalorimetry  analysis 
using  sputter  deposition.  The  thickness  of  each  layer  was  determined  to  ensure  a 
stoichiometric  composition:  the  samples  initially  had  a  3:1  molar  ratio  of  Zr:B4C  to 
produce  a  2:1  molar  ratio  of  ZrB2:ZrC  after  reaction.  Multilayer  samples  with  bilayer 
periods  of  5  nm,  10  nm,  and  a  total  thickness  of  80  nm  were  deposited. 

Fig.  10(a,b)  shows  the  enthalpy  rate  per  unit  area  as  a  function  of  temperature  for  ZrTEUC 
reactive  multilayers  with  bilayer  periods  of  (a)  5  nm  and  (b)  10  nm,  and  for  various 
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Fig.  10.  Enthalpy  rate  per  unit  area  vs.  temperature  for  Zr/B4C  multilayer  samples  with  a  bilayer 
period  of  (a)  5nm  and  (b)  10  nm.  (c)  XRD  results  on  the  as-deposited  multilayer  as  well  as  the  samples 
heated  to  1110  and  1500K.  (d)  TEM  micrographs  for  the  sample  heated  to  1500K:  (left)  a  bright  field 
image  with  FFT  analysis  for  selected  regions  and  (right)  diffraction  pattern. 


scanning  rates.  As  found  in  Zr/B  multilayers,  the  reaction  in  the  Zr/B4C  multilayers  also 
proceeds  in  two  stages  with  corresponding  exothermic  peaks:  the  step  that  corresponds  to 
the  first  exotherm  starts  around  400  K  for  all  samples  tested.  The  reaction  completes 
below  1000  K  for  the  samples  with  a  5  nm  bilayer  period  and  below  1150  K  for  the 
samples  with  a  10  nm  period.  Higher  heating  rates  and  bilayer  periods  result  in  higher 
peak  temperatures  for  the  first  reaction. 

The  second  stage  is  completed  at  approximately  1400  K  for  the  samples  with  a  5  nm 
bilayer  period  and  at  1550  K  for  the  10  nm  samples.  Again,  higher  heating  rates  and 
larger  bilayer  periods  shift  the  peak  temperature  to  higher  values,  but  the  effect  is  less 
pronounced  than  for  the  first  stage. 

To  determine  the  evolution  of  the  multilayer  structure  during  the  two  exothermic 
processes,  XRD  was  performed  on  multilayer  samples  with  10  nm  bilayer  periods  in  the 
as-deposited  state,  heated  to  an  intermediate  temperature,  and  after  completion  of  the 
second  exothermic  reaction.  As  shown  in  Fig.  10(c),  the  as-deposited  sample  has  a 
diffraction  peak  corresponding  to  crystalline  Zr  with  a  (0002)  texture  (Powder  Diffraction 
File  #050665).  No  other  reflections  were  identified,  indicating  that  the  as-deposited 
multilayers  consist  of  crystalline  Zr  and  amorphous  B4C.  The  intermediate  sample  has  a 
broad  peak  corresponding  to  an  amorphous  structure,  i.e.,  the  first  exotherm  is  caused  by 
the  intermixing  and  amorphization  of  the  two  constituent  layers.  The  final  sample  has  a 
strong  peak  corresponding  to  the  (0001)  ZrB2  reflection  (Powder  Diffraction  File 
#050665),  and  a  very  small  peak  caused  by  the  (1010)  ZrB2  or  the  (111)  ZrC  (Powder 
Diffraction  File  #050665)  reflections  -  clearly,  the  second  stage  corresponds  to  the 
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crystallization  of  the  amorphous  phase  to  form  a  supersaturated,  crystalline  ZrB2  phase. 
TEM  performed  on  the  same  sample  confirms  the  XRD  findings  (Fig.  10(d)). 

The  reaction  enthalpy  can  be  calculated  by  integrating  the  heat  flow  curves  in  Fig.  10 
with  respect  to  time.  The  reaction  enthalpy  is  independent  of  the  heating  rate,  but  varies 
slightly  with  the  multilayer  period  because  of  intermixing  that  occurs  during  sample 
deposition  -  the  total  thickness  of  the  intermixed  layer  is  approximately  1.9  nm  per 
bilayer  period. 

The  diffusion  model  developed  in  the  previous  section  can  be  used  to  analyze  the  Zr/E^C 
inter-diffusion  process.  According  to  Fig.  11(a),  the  activation  energy  for  inter- diffusion 
is  0.37+  0.08  eV.  The  crystallization  kinetics  can  be  characterized  using  a  Kissinger 
analysis  of  the  crystallization  peak  temperature  (Fig.  11(b)),  yielding  an  activation  energy 
of  1.85  eV.  The  activation  energies  for  diffusion  and  crystallization  obtained  for  the 
Zr/E^C  multilayers  are  both  smaller  than  for  the  Zr/B  multilayers.  Since  the  multilayered 
structure  were  synthesized  in  high  vacuum,  the  result  clearly  suggests  that  carbon 
enhances  the  formation  of  ZrB2  not  only  by  removing  surface  oxides,  but  also  by 
lowering  the  energy  barriers  required  to  form  the  compound.  The  lower  energy  barriers 
may  be  attributed  to  the  larger  electronegativity  difference  between  Zr  and  C  (1.22)  than 
between  Zr  and  B  (0.71),  which  leads  to  a  larger  binding  energy  between  Zr  and  C  than 
between  Zr  and  B  -  C  atoms  can  pull  Zr  atoms  toward  the  amorphous  region  more  easily 
than  boron  leading  to  more  distortion  of  Zr  lattice.  This  process  enhances  amorphization 
of  the  Zr  lattice  and  makes  more  room  for  B  or  C  atoms  to  diffuse.  This  last  notion  is 
explored  in  more  detail  in  the  next  section. 

A3.  Ab-initio  calculations  for  solid  state  Zr-B  alloying 

In  the  previous  section,  we  have  shown  that  inter-diffusion  of  B  and  Zr  is  facilitated  by 
the  formation  of  an  intermediate  amorphous  phase.  To  achieve  a  more  fundamental 
understanding  of  inter-diffusion  of  Zr-B,  we  conducted  a  first-principles  investigation  of 
the  inter-diffusion  process. 


(a)  (b) 


Fig.  11.  Kinetic  analysis  on  the  reaction  in  Zr/B4C  multilayer:  (a)  In ((//  +  H0)dH fdt)  versus  1  /kBT  for 
the  temperature  range  corresponding  to  interdiffusion  between  Zr  and  B.  Slopes  in  this  plot  yield  the 
activation  energy  of  the  diffusion,  (b)  Kissinger  analysis  for  crystallization  peaks  from  the 
measurements  on  the  samples  with  the  5  nm  period. 
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Fig.  12.  (a)  Positions  and  local  environments  of  the  interstitial  sites  in  Zr.  (b)  Activation  energies  of  a  single 
B  atom  along  different  diffusion  trajectories  calculated  by  the  NEB  method.  The  energy  profile  at  the 
starting  point  takes  into  account  the  probability  of  B  being  filled  at  the  initial  position  of  the  each  path. 


We  first  study  the  energetics  of  B  insertion  into  crystalline  Zr  by  determining  the  stable 
positions  of  a  single  B  atom  in  hexagonal  Zr  (P63/mmc,  a  =  b  =  3.23  A,  and  c  = 
5.17  A35,36).  A  single  B  atom  is  inserted  into  a  super-cell  at  nonequivalent  sites  as  shown 
in  Fig.  12(a).  Sites  include  the  octahedral  site  (Oct),  the  crowdion  site  (C),  the  basal  bond 
center  (Be),  the  hexagonal  site  (Hex),  the  plane  center  (Pc),  the  tetrahedral  site  (Tet),  and 
the  substitutional  site.  To  compare  the  thermodynamic  stability  of  B  at  each  interstitial 
site,  we  calculate  the  formation  energy  per  B  atom  Emin)  using  Eq.  (5.12) 

Emin)  =  LEnB-zr  -  (54-m)  EZr  -  nEB\/n.  (5.12) 

Here,  we  take  the  energies  of  a  Zr  atom  ( EZr )  and  a  B  atom  (EB)  in  bulk  form  as  reference 
energies,  with  £nB-zr  being  the  total  energy  of  the  system  containing  n  B  atoms  in  a  cell 
that  contains  (54  -  m )  Zr  atoms  (m= 0,  except  when  B  is  on  the  substitutional  site,  in 
which  case  m  =  1).  Table  2  lists  the  calculated  values  of  the  formation  energy  at  each 
interstitial  site  along  with  associated  structural  parameters,  where  sv  is  the  volumetric 
strain  and  CN  represents  the  coordination  number.  The  formation  energies  at  the  different 
positions  show  that  B  insertion  at  the  Oct  site  results  in  the  most  thermodynamically 


Table  2.  Structural  and  energetic  features  of  a  single  B  atom  at  various  sites  in  the  Zr  lattice 


Boron  sites 

Position 

Effl  (eV) 

CN 

BV  (%) 

Stability 

Octahedral  (Oct) 

(0.33,  0.67,  0.25) 

-0.53 

6 

0.732 

Stable 

Crowdion  (C) 

(0.83,0.17,0.25) 

0.59 

6 

0.435 

Metastable 

Basal  bond  center  (Be) 

(0.67,  0.83,  0.50) 

0.82 

6 

0.840 

Metastable 

Hexahedral  (Hex) 

(0.00,  0.00,  0.50) 

1.45 

5 

1.021 

Metastable 

Plane  Center  (Pc) 

(0.33,  0.67,  0.50) 

2.18 

3 

0.807 

Metastable 

Tetrahedral  (Tet) 

(0.67,  0.33,  017) 

2.03 

4 

1.134 

Unstable 

Substitutional  (Sub) 

(0,  0,  0) 

1.86 

12 

-1.605 

Metastable 
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stable  structure.  The  C,  Be,  Hex,  and  Pc  sites  are  defined  as  metastable  positions  because 
the  B  atom  would  be  stable  at  the  initial  insertion  site  upon  relaxation,  but  the  formation 
energies  are  positive.  The  Tet  site  is  unstable  in  that  a  B  atom  placed  at  this  site  will 
move  to  the  nearby  hexahedral  site  during  relaxation.  As  a  reference,  we  also  calculate 
the  formation  energy  for  substitutional  B  in  Zr.  The  large  formation  energy  indicates  that 
substitution  of  Zr  by  B  is  energetically  unfavorable  and  may  not  occur  at  low 
temperatures.  Note  that  B  insertion  at  interstitial  sites  with  lower  formation  energies 
results  in  relatively  small  volumetric  strains  (ev)  and  large  coordination  numbers 
representing  a  less  repulsive  and  more  cooperative  local  environment. 

In  order  to  determine  the  diffusion  pathway  of  B  in  crystalline  Zr,  we  have  calculated  the 
activation  energies  along  various  diffusion  trajectories  using  the  nudged  elastic  band 
(NEB)  technique.  In  the  calculations,  all  atoms  in  the  initial  and  final  states  are  first 
relaxed  in  a  super-cell  with  fixed  volume  and  shape.  We  then  sample  the  diffusion 
pathway  consisting  of  three  intermediate  configurations.  The  NEB  path  is  constructed  by 
linear  interpolation  of  the  atomic  coordinates,  and  then  relaxed.  The  activation  energy 
versus  the  reaction  coordinate  is  illustrated  in  Fig.  12(b)  for  several  diffusion  paths. 
Table  3  lists  the  calculated  energy  barrier  (AEneb)  and  the  energy  corresponding  to  the 
occupancy  probability  of  B  at  the  initial  position  of  the  diffusion  path  (AEj).  When  the 
initial  position  is  Oct,  A Ei  is  set  to  zero.  The  actual  diffusion  barrier  A Et  is  then 
AE{  +  AEneb.  The  NEB  simulations  indicate  that  the  Oct-C-Oct  transition  has  the  smallest 
diffusion  energy  barrier  (1.39  eV).  The  Oct-Bc-Oct  transition  is  the  next  preferred  path 
with  an  energy  barrier  of  1 .64  eV.  Although  the  difference  in  the  energy  barrier  between 
the  two  favorable  paths  is  not  large  (0.25  eV),  the  corresponding  Boltzmann  factor,  exp(- 
0.25/kgE),  where  ks  is  the  Boltzmann  constant  and  T  is  the  temperature,  indicates  that 
diffusion  along  the  second  favorable  path  is  approximately  18  times  less  likely  than  along 
the  first  path  at  1000  K,  the  temperature  where  the  intermixing  between  B  and  Zr  actively 
occurs27.  Other  transition  paths  do  not  contribute  to  the  diffusion  kinetics  much  as  they 
have  higher  energy  barriers.  Clearly,  B  transport  in  crystalline  Zr  is  highly  anisotropic 
with  the  hop  from  the  octahedral  to  the  crowdion  site  dominating  the  diffusion  kinetics  in 
the  basal  plane. 


Table  3.  Activation  energies  calculated  by  the  NEB  method  for  a  single  B  atom  along  various 
diffusion  trajectories. 


Diffusion  paths 

AEneb  (eV) 

AE.(eV) 

AEt(eV) 

Oct-C 

1.39 

0 

1.39 

Oct-Bc 

1.64 

0 

1.64 

Oct-Hex 

2.07 

0 

2.07 

Oct  -  Pc 

2.80 

0 

2.8 

Be -Hex 

0.46 

1.6 

2.06 

C-Hex 

0.73 

1.39 

2.12 
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Fig.  13.  Structural  and  energetic  features  of  ZrB0.25<x<2-  (a)  Representative  geometries  of  the  ZrB* 
structures,  (b)  Enthalpy  of  formation  of  ZrB*  in  this  work  (circles)  and  the  values  of  crystalline 
references  (squares),  (c)  Volumetric  strain  (d)  Coordination  numbers  of  Zr-Zr,  B-Zr,  Zr-B,  and  B-B 
pairs  in  the  ZrB*  structures,  (e)  B  diffusivity  vs.  composition  at  temperatures  of  1000,  1500,  and 
2000  K  obtained  from  ab  initio  MD  simulations,  (f)  Red  circles  show  the  diffusion  activation  energies 
calculated  from  MD  simulations.  The  error  bars  are  calculated  by  linear  fitting  the  MSD  over 
different  time  ranges. 


We  now  investigate  the  energetics,  structural  features,  and  diffusion  kinetics  in 
intermediate  compounds,  ZrBx,  0.25  <v<  2.  We  incrementally  increase  the  B 
concentration  in  the  Zr  host  by  inserting  B  into  the  largest  voids  of  the  super-cell.  In 
crystalline  Zr,  the  centers  of  the  largest  voids  correspond  to  the  octahedral  sites.  In  Fig. 
13(a),  we  show  four  representative  structures,  ZrB05,  ZrB,  ZrB  1.5,  and  ZrB2.  It  is  clear 
that  B  insertion  disorders  the  Zr  lattice  drastically. 

For  each  configuration,  we  calculate  the  enthalpy  of  formation  by  using  Ef=  [Zwr-nB  - 
mEzr  -  nE^\l{m+n),  where  m  and  n  are  the  number  of  Zr  and  B  atoms,  respectively.  The 
dependence  of  Efon  composition  is  shown  in  Fig.  13(b).  Note  that  the  Oct  site,  the  only 
stable  site  with  a  negative  value  of  the  formation  energy  shown  in  Table  2,  has  an 
equivalent  number  of  Zr  host  atoms.  Thus  the  ZrB  phase  represents  the  case  where  all 
Oct  sites  are  filled  by  B.  This  structure  is  the  most  energetically  favorable  among  the 
metastable  compounds.  We  include  for  comparison  reference  data  points  that  correspond 

on  00 

to  experimentally  found  or  theoretically  predicted  Zr-B  crystalline  phases  ’  .  As 
expected,  most  of  the  crystalline  phases  have  a  smaller  Ef  compared  to  the  metastable 
structures  at  the  same  concentration.  One  exception  is  that  the  ZrB  phase  in  this  study  has 
a  smaller  formation  enthalpy  than  the  rock-salt  phase  reported  in  experiments. 

Fig.  13(c)  shows  the  evolution  of  the  volumetric  strain  as  a  function  of  the  B 
concentration.  The  slope  of  the  volumetric  strain  remains  very  nearly  constant  throughout 
the  entire  range.  We  fit  the  dependence  of  the  volumetric  strain  on  the  B  concentration 
using  a  linear  function,  and  obtain  the  partial  molar  volume  of  B  in  Zr,  CIb-zt  =  0.51xl0'29 
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m3.  The  large  volumetric  deformation  during  Zr-B  alloying  is  accommodated  by  the 
amorphization  and  the  gradual  rearrangement  of  the  structure  through  Zr-Zr  bond 
breaking  and  Zr-B  cluster  formation.  To  examine  the  bonding  mechanism,  we  measure 
the  evolution  of  the  atomic  coordination  at  various  concentrations.  A  physically 
meaningful  measure  of  atomic  coordination  is  needed  to  define  the  atomic  bonds.  We 
define  Zr-Zr  atoms  to  be  bonded  if  their  distance  is  within  1 10%  of  the  Zr-Zr  bond  length 
in  bulk  crystalline  ZrB2,  which  corresponds  to  a  largest  bond  distance  of  3.50  A. 
Likewise,  we  define  that  a  Zr  atom  is  bonded  with  a  B  atom  if  their  distance  is  smaller 
than  2.81  A,  and  that  two  B  atoms  are  bonded  to  each  other  if  their  bond  length  is  less 
than  2.02  A.  Fig.  13(d)  shows  the  average  values  of  the  CN.  The  Zr-Zr  coordination 
gradually  decreases  indicating  a  constant  breaking  of  the  bonds  between  the  host  atoms. 
Interestingly,  the  rates  of  bond  breaking  and  formation  change  discontinuously  at  the  ZrB 
composition,  indicative  of  the  more  drastic  structural  change  that  occurs  beyond  this 
composition.  More  specifically,  the  CN  for  the  B-Zr  pairs  stays  nearly  6  for  ZrBx<i, 
which  is  the  same  as  the  CN  for  the  Oct  site  (Table  2).  When  the  B  concentration 
increases  further,  the  CN  starts  to  decrease  indicating  that  metastable  sites  are  now  being 
filled.  We  find  that  most  of  the  additional  B  atoms  in  ZrBx>i  occupy  Hex-like  sites,  which 
facilitates  amorphization  of  the  structure. 

We  quantify  the  diffusivity  of  B  in  various  Zr-B  compounds  by  performing  ab-initio 
molecular  dynamics  (AIMD)  simulations  based  on  density  functional  theory  using  the 
eight  relaxed  structures  with  various  concentrations  (ZrBo.25<x<2).  AIMD  calculations  are 
performed  at  temperatures  of  1500,  1625,  1750,  1875,  2000,  and  2500  K.  The  super-cells 
are  allowed  to  dynamically  equilibrate  at  each  temperature  for  10,000  time  steps,  each 
step  corresponding  to  2  fs.  The  diffusivity  of  B  can  be  determined  from  the  mean  square 
displacements  of  the  B  atoms  using  Einstein’s  relationship.  The  diffusivity  is  determined 
at  several  different  temperatures  and  is  fitted  using  an  Arrhenius  equation.  Fig.  13  (e) 
depicts  the  B  diffusivity  for  several  B  concentrations  at  temperatures  of  1000,  1500,  and 
2000  K.  For  all  temperatures,  the  diffusivity  is  minimum  at  ZrB.  The  effect  of  the  B 
concentration  on  the  diffusivity  is  striking,  especially  at  the  lowest  temperature  (1000  K), 
with  a  difference  of  four  orders  of  magnitude  between  extreme  values.  This  is  a 
consequence  of  the  large  activation  energy  for  diffusion  in  the  stoichiometric  ZrB 
structure  (Fig.  13  (f))  -  B  fills  all  Oct  sites  in  the  structure  making  B  diffusion  through 
the  structure  difficult. 

With  the  formation  of  an  amorphous  structure  observed  in  both  the  reactive  Zr/B 
multilayer  samples  and  the  AIMD  simulations,  it  is  possible  to  directly  compare 
experiments  and  simulations.  The  diffusion  activation  energy  obtained  from  nano¬ 
calorimetry  measurement  on  Zr/B  multilayers  is  shown  in  Fig.  13(f),  along  with  the 
simulation  results.  The  AIMD  result  at  the  same  B  concentration  is  in  good  agreement 
with  the  experimental  value. 

B.  Oxidation  of  ZrB2-based  coatings 

Oxidation  reactions  are  often  characterized  using  thermogravimetric  analysis.  In  this 
technique  changes  in  the  mass  of  a  solid  are  measured  as  a  function  of  temperature  and 
time,  and  related  to  the  extent  of  reaction.  The  technique  works  well  for  bulk  materials, 
but  is  less  effective  for  coatings  because  these  need  to  be  removed  from  their  substrate  to 
avoid  the  dilution  effect  caused  by  the  substrate39,40.  This  extra  step  in  sample  preparation 
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can  be  problematic  and  is  sometimes  not  feasible  when  the  coating  is  very  thin  .  Instead, 
AC  nano-calorimetry  can  be  used  to  characterize  the  oxidation  behavior  of  very  thin 
coatings  by  measuring  the  change  in  heat  capacity  of  a  sample  during  the  reaction41.  We 
first  demonstrate  this  technique  by  evaluating  the  oxidation  kinetics  of  sputter-deposited 
thin  films  of  zirconium  in  air.  We  then  apply  the  same  technique  to  study  the  oxidation 
kinetics  of  sputter- deposited  thin-film  ZrB2,  Zr(B,C)2,  and  Zr-B-Nb  samples. 


Bl.  Oxidation  kinetics  of  Zr 

The  oxidation  kinetics  of  Zr  can  be  evaluated  by  measuring  the  evolution  of  the  sample 
heat  capacity  either  isothermally  or  as  a  function  of  temperature.  Fig.  14  (a)  depicts  the 
results  for  a  set  of  isothermal  measurements  performed  in  air  on  200  nm  zirconium  films. 
The  change  in  heat  capacity  can  be  attributed  solely  to  the  oxidation  of  zirconium,  since 
no  phase  transformations  occur  in  zirconium  or  Zr02  in  this  temperature  range. 

If  the  rate  of  oxidation  is  limited  by  diffusion  of  oxygen  through  the  oxide  layer,  growth 
of  the  oxide  layer  is  described  by  the  following  equation42, 


dAh  _  K 
dt  ~  Ah+hi 


(5.13) 


where  A h  is  the  thickness  of  the  oxide  scale,  K  the  parabolic  rate  constant,  and  ht  the 
oxide  thickness  at  the  onset  of  the  isothermal  segment  of  the  experiment.  The  evolution 
of  AC  is  then  described  by  a  simple  quadratic  equation, 

AC2  +  2ACC(  =  2^t,  (5.14) 

where  2  is  a  material  constant  given  by  2  =Ah/AC  and  Cj=hJX.  The  parameters  Q  and  K 
are  readily  determined  from  a  least-squares  fit  of  Eq.  (5.14)  to  experimental  A C-t  data. 
The  fitting  results  for  the  oxidation  of  zirconium  are  shown  as  solid  lines  in  Fig  14(a). 
Equation  (5.14)  agrees  well  with  the  experimental  heat  capacity  curves  up  to 
approximately  80%  of  their  maximum  value,  demonstrating  that  the  reaction  indeed 
follows  a  parabolic  rate  law  and  oxidation  is  diffusion-limited.  As  most  of  the  zirconium 


Fig.  14.  (a)  Change  in  heat  capacity  as  a  function  of  time  for  five  different  isothermal  oxidation 
experiments  on  Zr;  the  solid  lines  represent  fits  of  Eq.  (4.19)  to  the  experimental  data,  (b)  Logarithm  of 
the  parabolic  constant  K  as  a  function  of  reciprocal  temperature  for  the  isothermal  (triangles)  and 
scanning  (squares)  measurements,  along  with  linear  least  squares  fits  and  data  from  references. 
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is  consumed,  however,  the  oxidation  rate  decreases  below  the  parabolic  rate.  This 
decrease  in  oxidation  rate  is  caused  by  the  finite  thickness  of  the  zirconium  sample  and 
the  non-uniform  temperature  distribution/oxidation  rate  across  the  sensors4,5,20.  An 
Arrhenius  graph  of  K  is  shown  in  Fig.  14(b)  for  all  isothermal  measurements,  together 
with  values  reported  for  bulk  zirconium43,44.  The  value  of  the  rate  constant  obtained  in  the 
calorimetry  measurements  is  roughly  one  to  two  orders  of  magnitude  larger  than  the 
literature  values,  because  the  very  small  grain  size  of  the  samples  (~20  nm)  promotes 
grain  boundary  diffusion  through  the  oxide45.  A  linear  least  square  fit  of  the  data  yields 
an  activation  energy  of  £fl=0.60  ±  0.09  eV,  comparable  to  the  activation  energy  for  low- 
temperature  oxidation  of  bulk  zirconium  in  the  literature43. 

If  the  reaction  mechanism  is  known  a  priori  (e.g.,  diffusion-limited),  it  is  possible  to 
quantify  the  reaction  kinetics  by  performing  scanning  calorimetry  measurements  on  a 
single  sample,  instead  of  isothermal  measurements  on  several  different  samples46,47.  This 
approach  is  demonstrated  in  Fig.  15  (a),  which  shows  three  different  temperature  scans 
for  a  single  zirconium  thin  film:  Scan  1  was  performed  in  vacuum,  while  scans  2  and  3 
were  performed  in  air.  Scan  1  provides  the  heat  capacity  of  the  unreacted  zirconium 
sample  as  a  function  of  temperature;  scan  2  shows  how  the  heat  capacity  of  the  sample 
changes  as  it  oxidizes  during  the  temperature  scan;  scan  3  depicts  the  heat  capacity  of  the 
fully  oxidized  sample.  Scan  2  overlaps  with  the  scan  1  below  500  K,  then  rises  above 
scan  1  to  join  scan  3  at  810  K.  The  overlap  between  scans  1  and  2  at  low  temperature 
indicates  that  the  zirconium  film  is  well  protected  by  its  native  oxide  and  that  there  is  no 
measurable  oxide  growth  in  this  temperature  range.  The  overlap  between  scans  2  and  3  at 
high  temperature  indicates  that  the  sample  is  fully  oxidized  at  the  end  of  the  second  scan. 
The  change  in  oxide  thickness  during  scan  2  is  then  given  by 

Ah  =  Cf^(hox-hi).  (5.15) 

C3_C1 

In  this  expression,  hi  and  hox  are  the  initial  and  final  oxide  thicknesses,  while  Ci,  C2  and 
C3  represent  the  experimental  heat  capacities  obtained  from  the  respective  scans.  The 


Fig.  15  (a)  Results  of  scanning  measurements  on  a  Zr  film  in  vacuum  (scan  1)  and  air  (scans  2  and  3). 
The  heat  capacity  of  the  addendum  from  the  calorimeter  is  included  in  all  three  measurements,  (b) 
Change  in  zirconium  oxide  thickness  as  a  function  of  temperature. 
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thickness  of  the  oxide  scale  is  plotted  as  a  function  of  temperature  in  Fig.  15  (b).  If  the 
initial  oxide  thickness  is  negligible,  the  parabolic  rate  coefficient  K  can  be  calculated 
directly  from  Eq.  (5.13);  if  the  initial  oxide  thickness  is  significant,  independent 
measurement  of  hi  is  needed.  Fig.  14(b)  shows  an  Arrhenius  graph  of  K  obtained  from  the 
scanning  measurement  on  zirconium,  along  with  the  results  from  the  isothermal 
measurements.  The  figure  demonstrates  a  linear  relationship  over  the  temperature  range 
from  550  K  to  800  K.  The  scanning  measurements  are  in  good  agreement  with  the 
isothermal  measurements  and  yield  an  activation  energy  of  Ea  =  0.59  +  0.03  eV, 
validating  the  scanning  AC  nano-calorimetry  approach  to  evaluate  reaction  kinetics  using 
a  single  sample.  While  isothermal  and  scanning  measurements  provide  similar  results,  the 
latter  is  performed  on  a  single  sample,  thus  providing  an  efficient  method  of  evaluating 
reaction  kinetics. 

B2.  Oxidation  kinetics  of  ZrB2 

The  oxidation  rate  of  sputter-deposited  ZrB2  to  form  Zr02  and  B2O3  was  measured 
isothermally  at  several  different  temperatures  to  determine  the  reaction  rate  law.  Fig.  16 
shows  a  summary  of  the  results  for  200  nm  and  270  nm  samples.  The  samples  are  fully 
oxidized  and  the  final  value  of  the  heat  capacity  indicates  that  there  is  no  significant 
evaporation  of  B2O3.  Compared  to  Zr  oxidation  (Fig.  14(a)),  three  observations  are 
noteworthy:  (1)  The  oxidation  rate  of  ZrB2  is  -  not  surprisingly  -  much  slower  than  for 
Zr;  (2)  the  ZrB2  oxidation  curve  has  a  significant  linear  segment,  followed  by  a  parabolic 
segment;  and  (3)  the  linear  segment  becomes  more  pronounced  with  increasing 
temperature.  The  last  two  observations  are  readily  explained  by  the  Deal-Grove  model 
for  oxidation.  According  to  this  model,  the  evolution  of  the  heat  capacity  of  a  thin-film 
sample  is  described  by,4 

AC2  +  —AC  =  KDt  (5.16) 


Fig.  16  Change  in  heat  capacity  as  a  function  of  time  for  isothermal  oxidation  experiments  on  (a)  200  nm 
ZrB2  samples  and  (b)  270  nm  ZrB2  samples;  the  solid  lines  represent  fits  of  Eq.  (5.16)  to  the  experimental 
data,  except  for  the  871  K  curve,  which  is  fitted  by  Eq.  (5.17).  The  sensor  for  the  960  K  isothermal 
measurement  broke  before  the  sample  was  fully  oxidized. 
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where  K&  and  Kr  are  proportional  to  the  diffusion  coefficient  and  the  rate  constant  of  the 
interfacial  reaction,  respectively.  Equation  (5.16)  agrees  well  with  the  experimental  heat 
capacity  curves  up  to  approximately  70%  of  their  maximum  value.  At  greater  values,  the 
oxidation  rate  drops  below  the  Deal-Grove  model  because  of  the  finite  thickness  of  the 
ZrB2  samples  and  the  non-uniform  oxidation  rate  across  the  sensors.  Thus  oxidation  of 
ZrB2  is  limited  by  the  interfacial  reaction  and  transport  of  reactive  species  through  the 
oxide  layer.  During  initial  oxidation,  AC  «  Kr>/Kr  and  the  quadratic  term  in  Eq.  (5.16)  is 
negligible;  oxidation  is  then  controlled  solely  by  the  interfacial  reaction  and  the  oxidation 
curve  is  linear  with  time.  As  the  oxide  layer  grows  in  thickness,  diffusion  through  the 
oxide  eventually  becomes  the  rate-controlling  process.  Thus,  oxidation  of  very  thin  ZrB2 
coatings  is  controlled  mainly  by  the  interfacial  reaction.  The  higher  the  temperature,  the 
thicker  the  coatings  for  which  this  is  the  case,  since  the  third  observation  implies  that  the 
ratio  Kd/Kr  increases  with  temperature. 

An  Arrhenius  graph  of  and  is  shown  in  Fig.  17.  Linear  least- square  fits  to  the  data 
yield  activation  energies  of  £r=0.71  eV  and  £b=2.35  eV  for  the  interfacial  reaction  and 
diffusion,  respectively.  The  activation  energy  for  diffusion  is  larger  than  for  the 
interfacial  reaction,  in  agreement  with  the  observation  that  Kv/Kr  is  an  increasing 
function  of  temperature.  The  activation  energy  is  similar  to  that  for  diffusion  of  oxygen  in 
bulk  ZrC>2,49  suggesting  that  diffusion  along  the  grain  boundaries  as  observed  for  the 
oxidation  of  Zr  is  impeded  in  this  case. 

At  temperatures  below  900  K,  the  Deal-Grove  model  does  not  provide  a  good  fit  to  the 
oxidation  data.  Instead,  ZrB2  oxidation  follows  a  logarithmic  rate  law, 

AC  =  A\n(Bt  +  1),  (5.17) 

where  A,  B  are  rate  constants.  The  logarithmic  rate  law  is  commonly  observed  for  very 
thin  oxide  layers,  but  has  not  yet  been  reported  for  ZrB2.  According  to  the  Mott-Cabrera 
model50,  this  behavior  is  observed  when  electrons  tunnel  through  a  very  thin  oxide  layer 
and  the  associated  electric  field  aids  transport  of  charged  species  through  the  oxide.  Since 
KdIKr  decreases  with  decreasing  temperature,  diffusion  of  the  reactive  species  through 
the  oxide  becomes  ever  more  difficult,  suggesting  that  the  Mott-Cabrera  model  should 
indeed  apply  at  low  temperatures. 


Fig.  17.  Logarithm  of  the  KR  and  Ku  as  a  function  of  reciprocal  temperature  for  the  isothermal  scans 
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Fig.  18.  (left)  Change  in  heat  capacity  as  a  function  of  time  for  ZrB2-ZrC  sample;  (right)  optical 
micrograph  of  a  ZrB2-ZrC  sample  after  oxidation. 


B3.  ZrB2-ZrC  oxidation  in  air 

It  has  been  shown  experimentally  that  during  oxidation  of  HfC  a  compact  layer  of 
HfC>2-xCy  is  formed  that  acts  as  a  diffusion  barrier  for  oxygen51,52.  Carbon  in  this  layer 
prevents  crystallization  of  the  layer  and  slows  down  diffusion  of  oxygen  significantly. 
These  findings  suggest  that  coatings  consisting  of  ZrB2  and  ZrC  may  in  fact  have  better 
oxidation  resistance  than  pure  ZrB2.  To  examine  this  hypothesis,  samples  were  fabricated 
by  co-sputtering  from  Zr  and  B4C  targets  in  a  3:1  molar  ratio,  resulting  in  amorphous 
2  ZrB2.ZrC  thin  films.  Fig.  18  shows  a  typical  result  obtained  for  isothermal  calorimetry 
at  940  K.  Unlike  the  oxidation  curves  for  Zr  or  ZrB2,  the  heat  capacity  shows  a  steep  an 
initial  drop,  followed  by  two  linear  segments  -  clearly  not  conventional  oxidation 
behavior.  Examination  of  the  oxidized  coating  using  optical  microscopy  reveals  that  the 
sample  experienced  gross  delamination  and  spalling  during  oxidation,  explaining  the 
drop  in  heat  capacity.  Delaminated  areas  are  circular  in  shape  suggesting  that 
delamination  is  caused  by  the  formation  of  small  blisters,  likely  as  a  result  of  CO  or  CO2 
formation  during  the  oxidation  process.  Evidently  thin  coatings  of  ZrB2/ZrC  are  not  good 
oxidation  barriers. 

B4.  Zr-B-Nb  oxidation  in  air 

Previous  studies  have  shown  that  doping  ZrB2  with  NbB2  slows  down  diffusion  of 
oxygen  at  temperatures  above  1600  K,  thus  enhancing  the  oxidation  resistance53,54.  Here 
we  evaluate  how  the  addition  of  NbB2  modifies  the  oxidation  behavior  of  ZrB2  thin  films 
at  temperatures  below  1000  K. 

Thin-film  samples  with  different  ZrB2:NbB2  ratios  were  synthesized  using  magnetron 
sputtering.  Fig.  19  shows  the  results  obtained  from  isothermal  calorimetry  measurement 
at  900  K  in  ambient  atmosphere.  Contrary  to  high-temperature  oxidation,  the  addition  of 
NbB2  enhances  the  rate  of  oxidation,  and  the  rate  of  oxidation  increases  monotonically 
with  Nb  content.  Furthermore,  the  oxidation  curves  for  the  doped  samples  are  more 
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Fig.  19.  Change  in  heat  capacity  as  a  function  of  time  for  ZrB2  samples  doped  with  different 
concentrations  of  NbB2  and  oxidized  at  900  K. 


parabolic  compared  to  un-doped  ZrB2,  indicating  that  the  interfacial  reaction  rates  are 
higher  and  that  oxidation  is  diffusion-controlled.  These  results  imply  that  the  intrinsic 
oxidation  rate  of  NbB2  is  higher  than  that  of  ZrB2. 

To  elucidate  the  temperature-dependence  of  the  oxidation  behavior,  isothermal 
calorimetry  measurements  were  performed  on  (Zro.86Nbo.i4)B2  at  several  different 
temperatures.  Fig.  20(a)  shows  the  oxidation  curves.  In  contrast  to  un-doped  ZrB2,  the 
curves  become  more  parabolic  with  increasing  temperature,  implying  that  the  ratio  K^/Kr 
decreases  with  temperature.  Fitting  the  oxidation  curves  with  Eq.  (5.16)  yields  values  for 


(a)  (b) 

Fig.  20.  (a)  Change  in  heat  capacity  as  a  function  of  time  for  (Zr0.86Nb0.i4)B2  during  isothermal 
oxidation  at  several  different  temperatures;  (b)  Arrhenius  graph  of  kinetics  parameters  for  the 
oxidation  of  (Zr0.86Nb0.i4)B2. 
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Kr  and  Kb.  An  Arrhenius  graph  of  the  results  is  shown  in  Fig.  20(b).  Linear  least-square 
fits  yield  activation  energies  of  £r=2.47  eV  and  Eb=2.  1  eV  for  the  interfacial  reaction 
and  diffusion,  respectively.  The  rate  constant  for  the  interfacial  reaction  is  much  larger 
than  for  un-doped  ZrB2  and  increases  much  faster  with  temperature,  i.e.,  the  addition  of 
NbB2  expedites  the  interfacial  reaction  significantly.  While  the  diffusion  rate  constant  is 
larger  than  for  un-doped  ZrB2  in  this  temperature  range,  the  activation  energy  is  smaller 
suggesting  that  Nb  indeed  slows  down  diffusion  at  more  elevated  temperatures,  in 
agreement  with  prior  observations53. 

C.  Mechanical  behavior  of  ZrB2  thin  films 

The  mechanical  behavior  of  bulk  ZrB2  depends  sensitively  on  the  densification  process 
used  to  synthesize  the  material  and  the  resulting  microstructure1,55'58.  While  there  are 
many  reports  on  the  mechanical  behavior  of  bulk  ZrB2,  the  behavior  of  ZrB2  thin  films  is 
relatively  less  explored59,60.  We  have  measured  the  stress-strain  curves  of  sputtered  ZrB2 
thin  films  between  room  temperature  and  1016  K.  ZrB2  thin-film  samples  with  a  900  nm 
thickness  were  deposited  onto  micro-machined  silicon  frames  with  integrated  tungsten 
micro-heaters  using  magnetron  sputtering  and  deformed  inside  a  scanning  electron 
microscope  (SEM)  with  the  custom-built  tensile  tester61,62,  described  in  an  earlier  section. 
Experiments  at  elevated  temperature  were  performed  by  allowing  the  ZrB2  tensile 
specimens  to  equilibrate  with  the  tungsten  micro-heaters  embedded  in  the  micro- 
machined  silicon  tabs. 

Figs.  21(a-d)  show  an  X-ray  diffraction  spectrum,  as  well  as  several  TEM  cross-sections 
for  a  typical  as-deposited  ZrB2  film.  The  TEM  micrographs  clearly  show  that  the 
sputtered  ZrB2  film  consists  of  two  distinct  layers:  an  amorphous  layer  (Fig.  21(d)  and  a 
crystalline  layer  with  a  fine  columnar  structure  (Fig.  21(c)).  The  transition  from 
amorphous  to  crystalline  is  caused  by  an  increase  in  the  substrate  temperature  during  the 
sputter  deposition  process. 

Fig.  21(e)  shows  the  stress-strain  curves  of  the  ZrB2  thin  films  for  several  temperatures. 
The  curves  are  linear  indicating  elastic  deformation  up  to  final  rupture  without  any  plastic 
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Fig.  21.  (a)  XRD  result  of  a  sputtered  ZrB2  film  in  the  as-deposited  state;  (b)  low-magnification  TEM  image 
of  the  cross-section  of  an  as-deposited  film;  (c)  high-magnification  TEM  image  of  the  cross-section  of  the  top 
layer  and  (d)  the  bottom  layer;  figure  insets  are  diffraction  patterns  confirming  a  crystalline  structure  for 
the  top  section  and  an  amorphous  structure  for  the  bottom  section;  (e)  stress-strain  curve  of  sputtered  ZrB2 
thin  films  tested  at  various  temperatures. 
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Table  4.  Summary  of  experimental  results  at  several  temperatures 


Temperature  (K) 

Elastic  modulus 
(GPa) 

Calculated  elastic  modulus 
(GPa) 

Tensile  strength 
(MPa) 

293 

295  ±  17 

532 

1300  +  79 

791 

250  +  7 

505 

1480  +  42 

1016 

170  +  5 

502 

950  +  27 

deformation.  Films  tested  at  room  temperature  show  an  exceptionally  high  fracture  stress 
of  1300  ±  79  MPa  (Table  4)  compared  with  reported  flexural  strengths  of  350  -  723  MPa 
for  bulk  ZrB21,56,57,63,64.  It  should  be  noted  that  uniaxial  tensile  testing  is  not  a  common 
method  for  evaluating  the  mechanical  behavior  of  brittle  ceramics  such  as  ZrB2,  because 
the  strength  of  a  brittle  material  when  loaded  in  tension  is  greatly  reduced  by  the  presence 
of  flaws  in  the  material.  Even  so,  sputter-deposited  ZrB2  thin  films  show  a  significantly 
higher  fracture  stress  in  tension.  This  observation  suggests  that  these  samples  have  a 
much  lower  flaw  density  than  bulk  samples  because  of  the  unique  sample  fabrication  and 
mechanical  testing  methods.  Table  4  lists  the  strength  as  a  function  of  temperature  -  the 
strength  increases  initially  with  temperature  and  still  reaches  950  MPa  at  1016  K.  The 
stiffness  of  the  films,  on  the  other  hand,  is  much  lower  than  for  bulk  materials  and 
decreases  significantly  with  increasing  temperature.  The  low  stiffness  of  the  films  is 
attributed  to  the  amorphous  region  as  well  as  the  high  density  of  grain  boundaries  in  the 
crystalline  region. 

The  fracture  toughness  of  the  ZrB2  films  was  measured  at  room  temperature  by 
introducing  pre-cracks  of  different  lengths  at  the  center  of  the  sample  gauge  section  using 
a  focused  ion  beam  and  loading  the  samples  in  tension  until  rupture  (Fig.  22(a,b)).  The 
stress-strain  curves  for  samples  with  various  pre-crack  lengths  are  shown  in  Fig.  22(c). 
According  to  linear  elastic  fracture  mechanics65,  the  plane-stress  stress  intensity  factor  K 
for  a  crack  in  the  center  of  a  thin  plate  subject  to  a  remote  tensile  stress  a  perpendicular  to 
the  crack  can  be  calculated  from 

*  =  (£).  (5.18) 

where  W  is  the  sample  width,  a  is  half  the  crack  length,  and  f(a/W)  is  a  function 
approximated  by 

f{$) = JlMl)  I1  -  °-°25©2 + °-o6(£)1-  <5-i9> 

The  crack  propagates  when  the  stress  intensity  factor  is  equal  to  the  fracture  toughness, 
Kc,  of  the  material.  Thus  the  fracture  stress  er/of  a  tensile  specimen  with  a  pre-crack  can 
be  written  as 


af  = 


a-1/2. 


(5.20) 
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SM,n  (Half  flaw  size  a)*,/2  (nn)‘,/2 

Figure  22.  (a)  SEM  image  of  a  freestanding  ZrB2  thin  film  with  a  pre-crack  of  size  2a  at  the  center  of  the  film, 
and  (b)  a  magnified  view  of  a  typical  pre-crack;  (c)  typical  stress-strain  curve  of  the  ZrB2  thin  film  tested  at 
room  temperature  with  different  crack  lengths;  (d)  fracture  stress  as  a  function  of  pre-crack  size. 

Fig.  22(d)  shows  the  fracture  stress  as  a  function  of  crack  length.  As  expected,  the 
fracture  stress  is  inversely  proportional  to  the  square  root  of  the  flaw  size.  By  fitting 
Eq.  (5.20)  to  the  experimental  data,  the  fracture  toughness  Kc  is  found  to  be 
2.57  ±  0.01  MPaVm,  which  is  similar  to  values  reported  for  bulk  ZrB2  ceramics56,64,66'68. 
It  should  be  noted  that  this  value  is  the  fracture  toughness  under  plane- stress  conditions, 
rather  than  the  usual  plane- strain  conditions,  which  are  more  appropriate  for  bulk 
materials.  For  brittle  materials  like  ZrB2,  however,  the  difference  between  both  is  small. 

VI.  Conclusions 

Scanning  AC  nanocalorimetry,  a  technique  developed  at  the  onset  of  the  project  to 
perform  calorimetry  measurements  at  elevated  temperatures,  has  made  it  possible  to 
make  nanocalorimetry  measurements  over  a  broad  range  of  temperatures  and  scanning 
rates.  The  dynamic  range  of  the  technique  makes  it  ideal  for  analyzing  the 
thermodynamics  and  kinetics  of  phase  transformations,  solid-gas  reactions,  and  solid- 
state  reactions  in  thin-film  samples. 

Scanning  AC  nanocalorimetry  analysis  of  the  solid-state  reaction  in  Zr/B  multilayer 
coatings  shows  that  the  formation  of  ZrB2  proceeds  in  two  distinct  steps:  an  inter- 
diffusion/amorphization  step  and  a  crystallization  step.  The  diffusion  step  has  a  small 
activation  energy  and  starts  at  a  surprisingly  low  temperatures.  The  formation  of  carbon- 
doped  ZrB2  from  Zr/B4C  reactive  multilayer  coatings  follows  a  similar  path,  but  the 
presence  of  carbon  increases  boron  mobility  and  lowers  the  activation  barrier  for  both 
diffusion  and  crystallization.  DFT  simulations  are  in  general  good  agreement  with  the 
experimental  findings. 
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The  oxidation  resistance  of  diboride-based  coatings  was  evaluated  by  measuring  the  heat 
capacity  of  a  sample  during  the  oxidation  reaction,  a  technique  first  demonstrated  on 
zirconium  thin  films.  Heat  capacity  measurements  were  performed  both  isothermally  and 
in  scanning  mode,  and  results  from  both  approaches  were  in  good  agreement.  Oxidation 
of  Zr  was  shown  to  be  diffusion-limited  with  an  energy  barrier  of  0.59  +  0.03  eV, 
suggesting  a  grain  boundary  diffusion  mechanism.  Oxidation  of  ZrB2  thin  films  followed 
the  Deal-Grove  oxidation  model,  with  activation  energies  of  0.71  eV  for  the  interfacial 
reaction,  and  2.35  eV  for  the  diffusion  step.  In  contrast  to  pure  ZrB2  films,  carbon-doped 
ZrB2  films  experienced  gross  delamination  and  spalling  as  a  result  of  gas  formation 
during  the  oxidation  process.  Doping  of  ZrB2  coatings  with  Nb  resulted  in  an  increase  in 
both  the  rates  of  the  interfacial  reaction  and  of  oxygen  diffusion,  lowering  the  overall 
oxidation  resistance  of  the  coatings.  However,  the  temperature-dependence  of  the 
diffusion  rate  constant  suggests  that  Nb  doping  may  be  beneficial  at  more  elevated 
temperatures. 

Tensile  experiments  performed  on  un-doped  ZrB2  thin  films  at  temperatures  ranging  from 
ambient  to  1016  K,  showed  an  ambient-temperature  tensile  strength  in  excess  of  1.2  GPa, 
far  exceeding  the  tensile  strength  of  bulk  ZrB2.  The  stiffness  of  the  films  was  lower  than 
for  bulk  materials  and  decreased  significantly  with  increasing  temperature.  The  fracture 
toughness  of  the  ZrB2  films  was  2.57  +  0.01  MPaVm,  which  is  similar  to  values  reported 
for  bulk  ZrB2. 
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VII.  Appendices 

A.  Thermal  analysis  for  nanocalorimetry  measurements  of  solid-state 

REACTIONS 

While  scanning  AC  nanocalorimetry  allows  direct  determination  of  the  heat  capacity  of  a 
sample,  it  is  also  possible  to  extract  information  on  the  enthalpy  rate  H  resulting  from  any 
solid-state  reaction  in  the  sample  by  considering  the  DC  component  of  the  calorimetry 
signal.  In  this  appendix,  we  describe  a  simple  but  effective  procedure  to  perform  this 
analysis.  Consider  the  following  energy  balance  equation, 

Po  =  Cp^  +  Lc(T0,x )  +  Lr(T0 )  +  H0(.T0,X )  (Al) 

Here,  Po  is  the  power  supplied  to  the  heater,  Cp  is  the  heat  capacity  of  the  control  volume 
(CV  -  i.e.,  sample  and  addendum).  Hq  and  T0  represent  the  reaction  enthalpy  rate  and 
temperature  of  the  CV  averaged  over  one  oscillating  period.  Lc  represents  the  power  lost 
to  the  environment  by  conduction  through  the  membrane  and  the  heating  element,  while 
Lr  represents  the  radiative  heat  loss,  both  quantities  appropriately  averaged  over  one 
oscillation  period.  Since  measurements  for  solid-state  reactions  are  typically  performed  in 
high  vacuum,  there  is  no  convective  heat  loss.  The  period- averaged  past  temperature 
history  of  the  CV  is  represented  by  The  analysis  procedure  is  based  on  the 
observations  that  the  radiative  heat  loss  is  a  function  of  temperature  only,  while  the 
conductive  heat  loss  depends  on  both  temperature  and  past  temperature  history. 
Furthermore,  for  small  temperature  oscillations,  the  averaged  heat  loss  terms  do  not 
depend  on  the  amplitude  of  the  temperature  oscillations.  The  reaction  enthalpy  rate,  H0, 
generally  depends  on  both  temperature  and  temperature  history. 

Consider  a  typical  experiment  that  consists  of  two  subsequent  calorimetry  measurements 
with  similar  but  not  identical  heating  rates.  The  sample  undergoes  an  irreversible  solid- 
state  reaction  during  the  first  measurement,  but  not  during  the  second.  The  energy 
balance  equations  for  both  scans  are  then 

pd1)(7o1))  =  CP1}  +  lc{t^.XW)  +  Lr(tV)  +  H0(rf\xw),  (A2) 

=  C®  ^  +  ^(r®,*®)  +  LR(rt ’®),  (A3) 

where  the  superscripts  (1)  and  (2)  refer  to  the  respective  measurements.  Taking  the 
difference  between  both  equations  for  T0  =  ,  and  rearranging  the  terms  leads 

to  the  following  expression  for  the  reaction  enthalpy  rate. 

H0(T0,xW)  =  (Po(1)<To)  -  P0(2)<To))  -  (c®  (^)(1)  -  C®  (^)(2))  - 

(A(W(1))  -  £c(W(2)))  =  AP  -  A  (Cp  -  A Lc,  (A4) 

in  which  the  radiative  heat  loss  term  has  been  eliminated.  Each  of  the  terms  in  the  RHS  of 
Eq.  (A4)  is  readily  evaluated  in  a  typical  calorimetry  experiment.  The  first  term 
represents  the  difference  in  power  supplied  in  the  two  scans  and  is  calculated  from  the 
applied  current  and  resistance  of  the  heating  element.  The  second  term  arises  because  the 
enthalpy  of  the  samples  accrues  at  different  rates  in  the  two  scans.  Since  scanning  AC 
nanocalorimetry  measurements  provide  the  heat  capacity  of  the  CV  as  a  function  of 
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temperature,  this  term  is  also  easily  calculated.  Evaluation  of  this  term  represents  a  clear 
advantage  of  scanning  AC  nanocalorimetry  over  straight  DC  nanocalorimetry  where  this 
term  is  much  more  difficult  to  determine  unless  the  two  scans  are  performed  at  exactly 
the  same  heating  rate  in  a  power-compensated  scheme  and  the  heat  capacity  of  the  CV 
does  not  change  after  the  reaction.  The  last  term  represents  the  difference  in  conductive 
heat  loss  as  a  result  of  the  different  temperature  histories  of  the  CV  during  the  two  scans. 
This  term  is  well  approximated  by  the  following  expression, 

Lc(T0(t),xw)  ~  Lc(T0(t),X™)  = 

-(2kmAm)  £  (£  (r0(V)  -  Km(t  -  t')dt')  \y=0  - 

(2fciA)£(Jo  (^oV)  -  r0®(t'))  Jfh(t  -  t')dt')  |*=o  .  (A5) 

where  the  subscript  m  refers  to  membrane  properties,  while  h  refers  to  properties  obtained 
by  appropriately  averaging  the  heater  and  membrane  properties.  A  detailed  derivation  of 
this  expression  and  the  definition  of  the  relevant  parameters  are  found  in  the  Appendix  B. 

In  a  typical  experiment  the  current  supplied  to  the  sensor  during  the  two  scans  is  ramped 
up  at  the  same  rate  and  the  difference  in  scan  rates  is  caused  by  the  enthalpy  associated 
with  a  reaction.  Consequently,  the  first  term  in  Eq.  (A4)  arises  because  the  resistance  of 
the  heating  element  is  a  function  of  temperature.  This  term  is  usually  the  dominant  term 
in  the  expression  for  the  reaction  enthalpy  rate.  The  difference  in  conduction  heat  loss 
arises  solely  as  a  result  of  the  slightly  different  thermal  histories  of  the  scans  and  is 
typically  quite  small;  the  same  is  true  for  the  term  associated  with  the  enthalpy.  The 
radiative  heat  loss,  which  at  elevated  temperature  is  often  the  largest  term  in  Eq.  (A2),  is 
automatically  eliminated,  although  some  care  needs  to  be  exercised  if  the  solid-state 
reaction  changes  the  emissivity  of  the  CV.  This  effect  shows  up  as  a  I4 -dependence  of 
the  reaction  enthalpy  rate  at  high  temperatures  and,  if  necessary,  can  be  eliminated  by 
fitting  the  experimental  data  over  a  temperature  range  where  the  enthalpy  rate  is  expected 
to  be  zero.  The  advantage  of  the  use  of  Eq.  (A4)  for  solid-state  reactions  then  lies  in  the 
fact  that  for  a  well-planned  experiment,  the  difference  terms  associated  with  the  heat 
capacity  and  conductive  heat  loss  are  small  compared  to  AP,  which  is  directly  determined 
from  the  sensor  output  signals. 


35 

Distribution  A  -  Approved  for  Public  Release 


B.  Thermal  analysis  for  conductive  heat  loss 

In  this  appendix,  we  provide  an  estimate  of  the  conduction  losses  through  the  membrane 
and  the  heating  element.  Let  y  be  the  coordinate  in  the  direction  perpendicular  to  the 
heating  element.  The  temperature  distribution  in  the  membrane  is  then  described  by  the 
one-dimensional  thermal  diffusion  equation4, 

pcPf=feg-2f(^-V),  (Bl) 

with  conditions 

T(y,  0)=  T0, 

T(0,t)=  T0+f(t),  (B2) 

T(co,t)  =  T0, 

where  r,  cp,  k,  e,  and  h  designate  the  mass  density,  heat  capacity,  thermal  conductivity, 
emissivity,  and  thickness  of  the  membrane;  a  is  the  Stefan-Boltzmann  constant.  The 
function  fit)  describes  the  temperature  history  of  the  heating  element.  Defining 


t  =  T-T0, 

(B3) 

k 

a  =  — , 

PCp 

(B4) 

p  _  8 £(JT03 

P  ~  hpcv  ’ 

(B5) 

and  performing  a  linear  expansion  of  the  radiation  term  in 
equation 

Eq.  (Bl)  leads  to  the  following 

dr  32t  n 

Tt  =  a^-PT- 

(B6) 

with  conditions 

r(y,  0)  =  0, 
t(0,  t)  =  f{t), 
t(oo,  t)  =  0. 

(B7) 

The  solution  of  this  equation  can  be  written  as4 

t(y,t)  =  fit')K(t  -  t')dt'. 

(B8) 

where 

K(t  t,,  _ 

y  J  2  Vi^(t-t')3/2 

(B9) 

The  integral  in  Eq.  (B8)  is  difficult  to  calculate  numerically  for  arbitrary  functions  f(t). 
Here  we  describe  a  simple  method  to  evaluate  the  integral  for  an  experimentally 
measured  temperature  history  based  on  the  realization  that  the  integral  is  analytical  for  a 
linear  temperature  history.  Specifically,  if  the  experimental  temperature  history  consists 
of  a  set  of  n  data  points  (4  Tj)  and  is  represented  by  a  linear  spline,  the  temperature 
profile  at  time  4  is  well  approximated  by 
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t (y,  t )  =  2?=12  j^i+104i  +  t')K(t  -  t')dt'  +  £jAn- !  +  Bn_±t')K(,t  -  (BIO) 

where 


and  Bi  = 


ti+i  ti  ti+1  ti 

The  heat  lost  by  conduction  into  the  membrane  is  then  given  by 
Lc  =  2kA(£?r?It  +  /„_!), 

where  A  is  the  total  cross-sectional  area  of  the  membrane  and 


(B 11) 

(B12) 


h  = 


+  Bjtn )  |~7===~  -  +  V/4"  (er/(7^(tn  -  t;+1))  -  erf{^{tn  -  t^))J  + 


Voir  l  y/tn~ti+ 1  Jt 

%  {erf(JP (fn  -  ti+i))  “  er/(V0(tn-t;))). 


2yfa(3 


(B13) 


-(^n-l  +  +  'fP*erfUP(tn  ~  tn-l))}  ~  ~  Ol-l))- 


(B14) 

Substituting  Eqs.  (B13)  an  (B14)  into  Eq.  (B12)  results  in  an  expression  that  can  be  used 
to  estimate  the  conductive  heat  loss  through  the  membrane  at  a  time  4  from  the  prior 
temperature  history  of  the  heating  element.  Evaluation  of  this  expression  requires 
knowledge  of  the  parameters  a  and  b,  which  are  readily  determined  experimentally4.  The 
main  approximation  in  the  analysis  is  the  linear  expansion  of  the  thermal  radiation  loss  in 
Eq.  (B6).  This  approximation  is  acceptable  as  long  as  the  neglected  radiation  loss  terms 
are  small  compared  to  the  other  terms  in  the  equation,  i.e.,  t  «  T(/bt,  where  t  is  the  time 
scale  of  the  experiment.  Given  the  approximate  character  of  the  analysis,  the  preferred 
method  of  eliminating  the  conduction  loss  from  experimental  nanocalorimetry  data  is  to 
take  the  difference  between  two  subsequent  scans  with  similar  thermal  histories.  The 
above  analysis  can  then  be  applied  as  a  small  correction  to  account  for  the  difference  in 
conduction  loss  as  a  result  of  the  slight  difference  in  thermal  history  between  the  two 
scans. 
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C.  ALPHA  -  ENTHALPY  ANALYSIS 

While  AC  nanocalorimetry  is  a  robust  method  for  measuring  the  heat  capacity  of  very 
small  samples  over  a  broad  range  of  heating  rates,  this  capability  comes  at  a  cost:  the  AC 
technique  does  not  provide  direct  information  on  the  enthalpy  released  during  a  phase 
transition  or  chemical  reaction  -  there  is  no  direct  correlation  between  the  peaks  in  the 
calorimetry  traces  and  the  enthalpy  of  the  corresponding  changes  in  the  sample.  Under 
some  limited  conditions,  however,  it  is  possible  to  relate  the  AC  signal  to  the  enthalpy 
production  during  a  phase  transition  or  reaction  using  the  analysis  described  below. 

Consider  the  case  where  a  reaction  or  phase  transition  occurs  during  the  calorimetry 
measurement.  For  a  single-step  process,  it  is  common  to  assume  that  the  rate  of  reaction 
or  transformation  follows  the  following  simple  relationship69 

^  =  fc(T)r(x)  =  J,  (Cl) 

where  a  represents  the  extent  of  conversion,  which  ranges  from  0  to  1,  and  where  Q  is  the 
total  enthalpy  of  reaction.  The  enthalpy  of  reaction  and  the  associated  rate  of  enthalpy 
production  H  are  taken  positive  for  an  endothermic  process.  The  extent  of  conversion  a  at 
time  t  can  be  obtained  as  the  fraction  of  the  total  enthalpy  of  reaction  released  at  that 
time69, 


x(t)  = 


Cm1 


(C2) 


Equations  (Cl)  and  (C2)  can  be  used  to  find  the  reaction  rate  and  the  extent  of  conversion 
from  the  rate  of  enthalpy  production  in  an  adiabatic  scanning  measurement. 

For  a  scanning  AC  measurement,  the  power  supplied  to  the  calorimeter  sensor  consists  of 
a  DC  component  and  an  AC  component.  The  resulting  temperature  response  follows  Eq. 
(3.2)  and  the  energy  balance  for  the  sensor  can  be  written  as20: 

P0  +  Pg  =  c  +  L(T0  +  e)  +  H(T0  +  e,  x ),  (C3) 


where  Po  and  Pe  represent  the  DC  and  AC  components  of  the  input  power,  respectively. 
If  6  is  small  and  the  angular  frequency  of  the  AC  component  is  large  enough  that  the 
change  in  a  during  a  single  oscillation  period  is  negligible  (k(T0)/a)  «  1),  H  can  be 
linearized  within  one  oscillation  period  so  that 


H  =  H0  +  di- 

u  dT 


To 


6  =  H0  +  aHd , 


(C4) 


where  H0  is  the  enthalpy  flow  for  the  reaction  at  temperature  To.  Substituting  Eq.  (C4) 
along  with  Eq.  (3.5)  into  Eq.  (C3),  and  separating  the  DC  and  AC  components  leads  to 
the  following  two  equations 


P0  =  Cp^  +  L0  +  H0,  (C5) 

Pe  =  Cp  —  +  ( oc0  +  aH)0.  (C6) 


If  the  heat  loss  during  the  measurement  is  negligibly  small  or  if  the  heat  loss  term  is 
eliminated  by  performing  differential  measurements70,71,  Eq.  (C5)  can  be  used  to 
determine  the  enthalpy  production  H0  during  a  reaction.  If  the  heat  loss  is  not  negligible, 
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Eq.  (C6)  can  be  used  to  obtain  information  on  the  enthalpy  production  instead.  Both  the 
heat  capacity  C  and  a  +  aH  can  be  determined  from  AC  calorimetry  data  using  the 
analysis  described  in  Section  B1  .In  Eq.  (C6),  aH  is  non-zero  only  during  the  reaction. 
As  a  result,  the  a  +  aH  curve  shows  a  well-defined  peak  during  the  reaction  and  aH  can 
be  obtained  from  the  a  +  aH  signal  by  baseline  subtraction.  The  quantity  aH  is  of  interest 
because  it  is  directly  related  to  the  kinetics  of  the  reaction.  Indeed,  combining  Eq.  (Cl) 
and  Eq.  (C4)  leads  to  the  following  expressions  for  a11 


a»  =  fr  =  <?£irrW- 

It  is  evident  from  this  equation  that 


aH 


a  In  k(T)/dT 


=  Qk(T)r(x)  =  H, 


(Cl) 

(C8) 


i.e.,  the  quantity  on  the  left  hand  side  of  this  equation  is  equal  to  the  enthalpy  production 
during  the  reaction.  Depending  on  the  precise  form  of  k(T),  this  expression  may  be  used 
in  Eq.  (C2)  to  determine  the  extent  of  reaction  during  an  AC  measurement. 

Equations  (Cl)  and  (C8)  can  be  used  to  obtain  information  on  the  kinetics  of  a  reaction  or 
transformation  from  a  series  of  scanning  DC  and  AC  nanocalorimetry  measurements.  For 
instance,  if  the  process  is  governed  by  Arrhenius  kinetics, 

k(T)  =  A  exp^-p  (C9) 

where  A  is  the  pre-factor,  Ea  is  the  activation  energy,  and  the  Boltzmann  constant,  then 
Eqs.  (Cl)  and  (C8)  can  be  rewritten  as 

\nH  =  lnAQ  +  \nr(x)-^~,  (CIO) 

ks  T 

In  aHT2  —  In  ^  =  In  AQ  +  In  (Cll) 

Equation  (CIO)  is  useful  in  the  analysis  of  adiabatic  DC  measurements,  while  Eq.  (Cll) 
can  be  applied  to  AC  measurements.  If  the  pre-factor  is  temperature  independent,  the  first 
term  on  the  right  hand  side  of  Eq.  (CIO)  is  constant.  Furthermore,  if  the  enthalpy 
production  is  evaluated  at  a  fixed  conversion  fraction,  the  second  term  is  constant  also, 
and  the  logarithm  of  the  enthalpy  production  changes  inversely  with  temperature.  Thus,  a 
graph  of  the  enthalpy  production  at  a  fixed  conversion  fraction  as  a  function  of  1  IT  yields 
a  straight  line  with  EJkB  as  slope.  Equation  (Cll)  is  a  special  case  of  Eq.  (C8)  and  shows 
that  anT2  is  proportional  to  the  enthalpy  production.  Thus,  the  conversion  fraction  for  the 
AC  measurement  can  be  calculated  from  Eq.  (C2),  where  olhT2  takes  the  place  of  H  as  the 
integrand.  A  graph  of  ln^z/T2)  as  a  function  of  1  IT  yields  a  straight  line  with  slope  EJkB, 
as  long  as  clmT2  is  evaluated  at  a  fixed  conversion  fraction.  The  right  hand  sides  of 
Eqs.  (CIO)  and  (Cll)  are  identical,  so  that  graphs  obtained  from  DC  and  AC 
measurements  are  offset  by  In (EJkB).  Since  the  activation  energy  is  not  known  a  priori, 
the  data  can  be  combined  into  a  single  graph  for  fitting  purposes  using  a  recursive 
approach.  The  combination  of  DC  and  AC  results  enables  kinetic  analysis  of  reactions 
over  a  wide  range  of  scanning  rates.  Neither  Eq.  (CIO)  nor  (Cll)  requires  a  constant 
scanning  rate,  but  rather  the  direct  measurements  of  H  and  aH.  This  is  an  especially 
useful  property  for  nano-calorimetric  measurements,  where  the  latent  heat  can  easily 
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change  the  scanning  rate,  and  feedback-control  of  the  temperature  is  not  always  feasible. 
The  same  technique  can  also  be  applied  to  non- Arrhenius  type  processes  as  illustrated  in 
the  following  Appendix. 

D.  Nucleation 

As  an  example  of  how  the  method  discussed  in  Appendix  C  can  be  used  to  investigate  a 
variety  of  phenomena,  we  describe  how  AC  and  DC  nanocalorimetry  can  be  combined  to 
evaluate  the  solidification  kinetics  of  thin  Bi  droplets  over  a  broad  range  of  cooling  rates. 

When  a  very  thin  coating  is  deposited  on  a  substrate  with  high  interfacial  energy  and  then 
melted,  the  coating  breaks  up  into  a  large  number  of  very  small  isolated  islands  to 
minimize  the  interfacial  energy  between  the  coating  material  and  the  substrate.  The 
coating  material  forms  a  dispersion  of  individual  droplets,  and  solidification  of  the 
dispersion  requires  nucleation  of  the  solid  phase  in  each  individual  droplet.  Under  these 
conditions,  solidification  is  nucleation-controlled  and  the  enthalpy  production  during 
solidification  is  given  by72 

H  =  Qk(T)r(x ),  (Dl) 

where 


k(T)  =  Aex  p(-jjjj),  (D2) 

r(x)  =  1  —  x.  (D3) 

In  these  equations,  which  derive  directly  from  classical  nucleation  theory73'75,  Q  is  the 
latent  heat,  k(T)  is  the  nucleation  frequency,  A  and  B  are  parameters  that  are 
approximately  constant74'77,  and  AT  is  the  degree  of  undercooling  of  the  islands  below 
their  melting  point.  Equation  (D2)  can  be  rewritten  as 

lnfc(T)  =  In  H-  ln[(l  -  x)Q]  =  \nA-^.  (D4) 

From  Eq.  (C8),  it  follows  that 

ln/cCT)  =  In  (|y  |)-ln[(l  -  x)Q]  -  In  B  =  (D5) 

where 


_  dH  _ 

aH  ~  dT  -  Q 

_  AT-2T 
9  T’2^7’3* 


dk(T) 

dT 


r(x) 


(D6) 

(D7) 


aH  is  directly  obtained  from  AC  calorimetry  results,  as  detailed  in  Section  Bl.  Equation 
(D4)  provides  an  expression  to  calculate  the  nucleation  frequency  k(T)  from  the  enthalpy 
production  in  an  adiabatic  DC  measurement,  while  Eq.  (D5)  can  be  used  to  compute  k(T) 
from  AC  measurements.  Fig.  (a)  shows  as  a  function  of  cooling  rate  the  undercooling 
obtained  for  a  200  nm  Bi  coating  evaporated  onto  the  silicon  nitride  surface  of  a 
calorimetry  sensor72.  The  undercooling  of  the  Bi  dispersion  changes  by  approximately 
20  K  over  three  orders  of  cooling  rate.  Fig.  (b)  shows  a  combined  graph  of  the  nucleation 
frequencies  for  DC  and  AC  cooling  scans  obtained  using  Eqs.  (D4)  and  (D5).  In  this 
graph,  the  value  of  B  needed  to  evaluate  the  nucleation  frequency  from  the  AC  data  was 
determined  using  recursive  least  squares  regression.  The  figure  clearly  illustrates  that  the 
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Fig.  (a)  Undercooling  of  Bi  at  50%  conversion  as  a  function  of  cooling  rate;  (b)  Graph  of  the  nucleation 
frequency  in  Bi  as  a  function  of  1/TAT2. 

AC  and  DC  data  are  in  good  agreement  with  each  other  and  that  the  slopes  of  both  sets  of 
data  are  the  same.  The  linear  relationship  between  the  data  indicates  that  a  single 
nucleation  mechanism  is  active  over  three  orders  of  magnitude  of  cooling  rate.  The 
kinetic  factor  B  and  the  pre-factor  A  can  be  calculated  from  the  slope  and  intercept  of  the 
linear  fit.  These  parameters  can  then  be  related  to  the  interfacial  energy  and  catalytic 
activity  of  the  interface  using  the  equations  of  classical  nucleation  theory. 
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